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ABSTRACT 


At the end of one year the program is on schedule. EFG ribbon growtii 
system JFL No. 2 is being used to perform meaningful growth rate and ribbon 
thickness experiments. The new, wide ribbon growth system is in operation. A 
theoretical study of stresses in ribbons is reported. EFG ribbons have been 
observed to exhibit a characteristic defect structure ^^bich is orientation dependent 
in the early stages of growth. 
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I. INTRODUCTION 


The objectives of this program are to produce at hi^ speed, thin, wide sheets 
of silicon of the quality needed to fabricate 10% efficient solar cells. This is part of 
a total Mobil Tyco Solar Energy Corporation objective to produce silicon ribbon for 
solar cells at a cost vdiich will allow their wide scale use for generation of electrical 
power. The edge-defined, film-fed growth process was developed for the commercial 
production of continuous shaped single crystals of sapphire from die melt and was 
applied to growth of silicon ribbon for solar cells partly under NSF Grant GI37067X 
via Harvard University, JPL Contract 953365, and under NSF Grant G143873. The 
basic feasibility of the application of EFG to the growth of silicon ribbon has been 
proven and the theoretical base for extending the technique to the efficient production 
of sheet silicon has been established. 

The two major tasks in this program to produce large area silicon shee by 
EFG are: the growth of thin ribbon at high speeds in existing apparatus an 1 the 
design and construction of a system for the growth of wide ribbons. These two tasks 
will be combined toward die end of the program to achieve the ultimate program goal. 
Concurrent with the development of the growth process and ribbon growth apparatus 
will be material and solar cell characterization programs and theoretical studies of 
various aspects of EFG silicon sheet growth including melt replenishment, the effect 
of h^^ and aspects of heat flow in the EFG silicon growth system. 

At the end of one year the program is on schedule. Work is progressing along 
sevoral fronts. Details of the progress achieved in the areas of thin/fast growth, 
"wide" system engineering and assembly, characterization, and theoretical studies 
follow in the Technical Discussion. 
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II. TECHNICAL DISCUSSION 

A. Apparatus 

One major task of this progrr m is to grow thin ^~0.1 mm) silicon ribbon at 
high growth rates ( 7.5 cm/min) . The?c goals are being pursued in a resistance 
heated EFG system designated JPL No. 2 (Fig. 1) . The system was modified at the 
beginning of this program to incorporate all current state-of-the-art EFG silicon 
ribbon technology. This modification included planning and building of a variable 
top heater assembly, the design of the power supplies and controller for the new top 
heater assembly, the design of specific crucible holder parts and an Improved 
crucible support pedestal assembly and a general debugging, rewiring and system 
alignment effort. 

The final configuration of the furnace chamber is virtually indistinguishable 
from the growth system previously described and illustrated,^^^ The puller is a 
ballscrew driven fixed stroke machine capable of pulling ribbon at up to 12.7 cm/min. 
The furnace itself is comprised of a vertical cylindrical graphite crucible heating 
element insulated with concentric carbon and graphite felt heat shields. Die heaters 
or 'trimmers* Include ribbon face heaters wired in parallel and end heaters con- 
trolled Independent of the face heaters and independent or one another. Molybdenum 
and graphite shields and the graphite fiberform package complete the top insulation. 
Real control of ribbon spreading and left-right temperature balance is achieved with 
the end heaters. All of this is packaged in a water-jacketed stainless steel furnace 
chamber. 

To achieve the thermal gradients necessary to grow low stress ribbons at 
high growth rates, an afterheater is currently being designed for this growth system. 
The layout is nearly complete. The heater will be flexible to provide various gradients 
as desired. 
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A new ribbon growth system, designated JPL No. 1 is design«jd to have the 
capacity to produce 125 cm long, 7.5 rm wide, 0.1 mm thick silicon ribbon at speeds 
of up to 7.5 cm/min. Experience with JPL No. 2 has contributed substantially to the 
design of JPL No. 1. We expect that new thin fast/growth technology developed on 
JPL No, 2 will contribute to experiment planning on the new mach ine as v/ell. 

JPL No. 1 was operational before the end of this rq>ortirg period (F't. 2). 

It is comprised of three elements; the puller, the furnace char.uer, and the 
electrical system. The puller is a 192 cm fixed stroke dovetailed slid ’Adth saddles 
which are ballscrew driven. The screw is driven by a closec' loop controlled torque 
motor. The furnace chamber is a horizontal, water-jacket-fxi, stainless steel (Fig, 3) 
can containing concentric semicylindrical molybdenum and tungsten heat shields (Fig. 
4), hi^ purity graphite heating element and crucible holdei*, and trough shaped quartz 
crucibles (Fig, 5a). The horizontal furnace geometry was selected so that the geom- 
etry Oi th.e physical system and of the thermal gradients more closely approximate the 
geometry of the product. The power is supplied to the main heating element from an 
SCR controlled 25 kW transformer. Independently controlled die heater /trimmers 
and an active afterheater are powered by an additional 10 kW supply (Fig. 5b). 

Though simple and dependable, the resulting system should be flexible enough 
to accommodate the host of experimental variables including 2.5 to 7.5 cm wide ribbon 
growth at speeds up to 7.5 cm/min as well as supplementary studies including contin- 
uous melt replenishment. 

The JPL No. 1 design and procurement efforts of the first three quarters of 
tills year have been reported. During the last quarter all the major components were 
received and ai^sembled. Details and debugging have been attended to. A seed holder 
was designed (Fig. 6a), radiation shields for the windows, water-cooling for top and 
bottom flanges and puller slide protection circuits have been designed and installed. 
SCR firing circuit and temperature controller malfunction have beai repaired. The 
main healing element has been heated up to 1800**C. Heat shield warping observed has 
be«i minimized by adding tungsten nuts and bolts and spacers. The die holder/ 
gradient controlling cartridge was assembled (Fig. 6b) . 

The system has been run up past 1400®C repeatedly. Temper aoires inside 
the shield package are essentially uniform. Excqjt for the ends of the molybdenum 
bolts, the outside of the shield package never exceeds ~ 850®C. A 300 g silicon 
charge was melted using less than 30% of the available 25 kW power. The third growth 
attempt yielded ribbon ~4.6 cm wide using less than 15 kW in the naln heater, die 
heaters, and afterheaters total. 
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Fig. 5. (a) Quartz 'trough shaped' crucible and graphite crucible holder. 

(b) JPL No. 1 control console containing closed loop motor speed, 
afterheater, and die heater controls and power supply and main 
heater cortrol. 
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Fig. 6. (a) x/y/0 pull head with microswitch push and pull protection, 
(b) Die holder-thermal trimmer cartridge. 
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B. Thermal Profiles and Thermal Stresses in the Growth of Rlhhon Crystals 
1. Introduction 

Attempts to grow ribbon-shaped crystals of silicon from the melt by 

techniques such as EFG or laser recrystallization have encountered problems with 

residual stresses in the grown crystals. Experience with EFG silicon ribbon growth 

has indicated that the stress problem became worse as faster growth rates and 

larger ribbon widths were sought. Since increasing the growth rate and ribbon width 

are of primary importance in developing the EFG process to attain its full potential, 

we have undertaken to consider the origin of the residual stresses in ribbon-shaped 

crystals and to identify the parameters that determine their magnitude. The analysis 

described here, of course, applies equally to other sheet growth techniques of silicon 

from the melt such as the Inverted Stepanov technique, the Ribbon-to-Ribbon process 

by laser recrystallization, and web-denontic growth; research on these techniques 

(as well as on EFG ribbon growth described here) is being currently sponsored by 

( 2 ) 

the Jet Propulsion Laboratory.' 

An identification of the origin of the residual stresses is made in Section 
it is indicated that the stresses result from the plastic deformation which the crystals 
undergo at the high temperatures. The plastic deformation, in turn, is caused by 
the thermal stresses which arise as a result of non-uniform temperature distributions 
in the ribbon. The problem of thermal profiles in ribbon-shaped crystals is con- 
sidered in Section 1I.B.3. The most useful approach to determining the temperature dis- 
tribution in the ribbon turns out to be numerical modelling on the computer; the model 
permits close interaction between the theoretical calculations and the design of 
experimental hardware. In Section 11.B.4, we apply thermoelasticity theory to predict 
the effects of the thermal profile on the thermal stresses in ribbon-shaped crystals. 

Of basic concern here are the two interrelated problems; one is to minimize the 
extent of plastic deformation and, therefore, the strain induced imperfections, the 
other is to reduce, to a tolerable level, the residual elastic stress in the crystal 
when it reaches room temperature. Various experimental approaches which are 
being undertaken to solve the thermal stress problem in EFG ribbon growth are dis- 
cussed in Section Il,B,4.d; a recent strategy which can lead to a minimization of both the 

plastic strain and the residual stress is also described. An oral version of the 

(3) 

following text was presented' at a seminar sponsored by the Jet Propulsion Labora- 
tory. 
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2. Origin of the residual stress 

It is worthwhile to review here briefly some basic concepts concerning 
residual stress in crystalline materials. Residual stress is a stress system satisfying 
internal equilibrium with no external loads or temperature gradients; it is accom- 
panied by an elastic strain distribution in the material. On a microscopic scale, 
residual stress is associated with the microstructure of the material such as dislo- • 
cations, inclusions, deformation twins, and dislocation pile-ups; these residual 
stresses are short range, but the defects play a fundamental role in plastic deforma- 
tion and in fracture of the material and in its electrical properties. Residual stresses 
on a macroscopic scale can be caused by mechanical effects (e.g., forced alignment of 
parts, or loads causing non-uniform plastic flow or creep), by thermal effects (e.g., 
thermal stress due to non-uniform changes in temperature causing plastic flow), and 
by chemical or compositional effects. The present analysis is concerned with the 
thermal stress in the growth of ribbon-shaped crystals, i.e., with die plastic strain 
induced imperfections and the residual stress in the grown crystal. Residual stresses 
can be generally relieved by small amounts of relaxation or plastic strain; if present, 
the stresses can cause or promote fracture, buckling, plastic flow and creep. Resi- 
dual stress in silicon can be measured by X-rays or birefringence, or, more typically, 
by observing the strain during removal of material (e.g., by scribing and splitting 
the ribbon crystals and measuring the curvature after splitting^^^ - see Fig. 7b). 

The residual stresses we have measuied in EFG silicon ribbons by such 
splitting technique were generally in the range 0 to 30, 000 psi (the fracture stress 
of silicon at room temperature is approximately 80, 000 psl^^^ ) . Cutting into solar 
cell blanks is generally difficult if the stress exceeds ~5, 000 psi; the uniformity 
of ribbon thickness also affects the ease of cutting. Spontaneous shattering of the 
ribbon during growth or handling, or excessive fracture during scribing is probably 
indicative of residual stresses in excess of 30, 000 psi. 

From the observed shapes of the split or buckled ribbons (see Fig. 7), we 

can already dcjuce some information on the elastic strain in the mateiial. In 

essence, the ribbon is longer in a central longitudinal section than along the edges; 

on a microscopic scale, the length variation across the ribbon width is the result 

of extra half-planes of dislocations. A simple calculation of the dislocation dr nsity 

for a typical stress of 15, 000 psi (which corresponds to a split of 0.52 in. over a 

4 

ribbon length of 20 in. in a 1 in. wld'^ ribbon) yields a density of 8 x 10 cm . In 
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Fig. 7 





(o) 



(b) 


. ( a) Portion of buckled ribbon showing ripples on ribbon surface, 
(b) Portion of scribed and split ribbon with a residual stress 
of ~7000 psi. 
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the analysis of the thermal stress in the ribbon, we will point out how the plastic 
deformation during growth, in response to the thermal stress, can lead to the in- 
crease in the lengtli at the ribbon center, 

Tlie basic question we must first consider is: vdiat is the origin of the resid- 
ual stresses in a non-uniform temperature gradient, and \^diat are die parameters 
that determine thdr magnitude. When the ribbon is cooled from the melting point 
to room temperature, each portion or element of the ribbon will contract propor- 
tional to its own temperature fall. If adjacent elements contract differently, they 
will restrain each other’s contraction, and hence stresses must result. External con- 
straints to thermal contraction can also cause thermal stress, but this does not 
apply to the crystal growth situations treated here. 

The effect of the temperature distribution along the ribbon length (or the 
"vertical" temperature profile) on the thermal stress is illustrated in Fig. 8. We 
consider the ribbon to consist of a number of narrower strips, as shown in Fig. 8a. 

If the ribbon is at a uniform temper? ture, the strips are parallel -sided (Fig. 8a) 
and there is no stress. If the ribbon is in a uniform temperature gradient, each 
strip tapers uniformly as a result of tliermal contraction (assuming, of course, that 
the thermal expansion coefficient is ind^endent of temperature). Figure 8b shows 
this, grossly exaggerated. Again, adjacent strips conform to each other without 
bending. Therefore, the ribbon is stress free, and will remain stress free as long 
as it is in the same temperature gradient. If the temperature gradient becomes non- 
uniform (i.e., non-linear), then the shape of each strip, if it were not attached to its 
neighbors, would be as shoiwi in Fig. 8c. It is apparent that adjacent strips no longer 
conform to each other without bending; the strips can be made to fit by lateral defor- 
mation in their own plane, thus thermal stresses must arise. It is also apparent that 
the stress will be proportional to the curvature or bending that must be imposed to 

fit the strips together. The curvature is, of course, proportional to the rate of change 

2 2 

of the temperature gradient, i.e., d T/dx , and also increases with increasing ribbon 
width.. The expllcii. dependence of the thermal stress components on the derivatives 
of the vertical temperature profile T(x) and on the riobon width will be derived in a 
later section. The important point to realize from this qualitative argument is that, 
in a realistic crystal growth situation, it is impossible to cool the ribbon to room 
temperature without developing stress in the ribbon at some point. 
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3. Thermal profiles in the growth of ribbon crystals 
a. Analytic calculation of thermal profile 

The first problem is to determine the temperature distribution in 
the ribbon under a variety of growth conditions. Extensive theoretical work was per- 
formed in a recently concluded NSF program^®^ to obtain analytical solutions for the 
vertical temperature profile and to relate the ribbon thickness to the gradioit in the 
meniscus and the growth rate. Numerical analyses of the horizontal (i.e., along the 
ribbon width) and transverse (i.e., throu^ the ribbon thickness) temperature profiles 
were also carried out. 

An analytic expression for the gradient in the vertical profile in the ribbon is 
17) 

given by the following: ' 

§■ ‘(Tf-r-V ■ [i + 41„(T/T„)]|' . (1) 

Here we have assumed that the heat flow is one-dimensional (i.e., along the vertical 
or x-direction) , the ribbon width is much greater than the ribbon thickness (t), the 
thermal conductivity is proportional to 1/T, and the radiating environment seen 
by the ribbon is at a uniform temperature Also, the heat carried by the moving 
ribbon has been neglected in calculating the steady-state temperature profile. The 
vertical temperature profile can be readily obtained by numerical integration of 
Eq. (l)withT(0) = T^. 

Equation (1 ), in conjunction with a steady- state heat balance condition at 
the solid-liquid interface, has been useful in relating the ribbon thickness to the 
growth rate and the gradient in the liquid; the effective environment temperature 
has been determined from correlations of experiment and theory. Otha: than 
describing the gradient at the interface, however, this analytic approach is insuffi- 
cient to describe the vertical profile in the ribbon to the extent that useful d^ermina- 


t The other terms in the above expression are the Stefan-Boltzmann constant cr, the 
emissivity of the ribbon surface c , and the thermal conductivity of solid silicon 
k at the melting temperature T . 
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tions of the thermal stresses can be made. The main problems ai*e the oversimplifi- 
cation of the radiating environment and the assumption of a stationary ribbon in the 
analytic calculations. The inclusion of a more realistic thermal environment (which 
can also be readily correlated to the crystal growth setup) and of the specific heat 
transport in the ribbon into the solution must perforce involve numerical analysis on 
the computer. The remainder of this section describes a computer model developed on 
a previous Mobil Tyco funded program, and examines the predictions of this model re- 
lating to the current program's objectives of fast, wide and stress -free ribbon growth. 

b. Numerical calculation of thermal profile 


Figure 9a depicts schematically the geometry used in the numerical 
calculation of the temperature profile. Again, only one-dimensional heat flow is con- 
sidered, i.e., the profile along the vertical (or crystal pulling) axis of the ribbon 
is calculated. The afterheater block dimensions Xg (hei^t above the interface) and 
yg (distance from ribbon surface), and the environment temperatures below, along 
and above the block are variable; a maximum of 10 temperature values can be speci- 
fied along the afterheater. "Tie steady-state temperature is calculated at 40 equally 
spaced incremental points along the ribbon; the length L, thickness and growth rate 
of the ribbon are variables in the program. 

The model considers conductive and convective heat transport in the ribbon, 
and radiative and convective cooling at the ribbon surface. The heat balance condi- 
tion for the ith element in the ribbon is given by (see Fig. 9b); 


Wt 

/T. , - 1 
k I 


^i \ Ax 

/ 1+1 \ Ax 

/J 

- 2Ax (w + 1) • 

^ RAD ^ CONV 
^i ^i 

= 


t. 


P CVg wt 



( 2 ) 


In the above expression, w is the ribbon width, C is the specific heat of solid silicon, 
Vg is the crystal growth rate, and k^ = 2 k^ ^ ^i-1^ average thermal 

conductivity for the i*^ and (i-1)^ elements. The radiative heat loss from the i*^ 
element (Q^^^^) is a function of T^ and € , and of the emissivities, temperatures 
and angle factors of the surfaces which comprise the environment (i.e., the after- 
heater block). The convective heat loss from the i^^ element is given by = 

h (T, - T ), where h = Nu • k /\q is the heat transfer coefficient, T„ is the effective 

* 8 th ^ ^ * 

gas temperature at the i element (assumed to be intermediate to T. and the after- 
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Fig. 9. (a) Schematic of thermal model used in the calculation of the 
vertical temperature profile in the ribbon, f W Heat flow 
terms at the i*^ element in the ribbon (qCOND - conductive heat, 
qSH s Mecific heat convection, qCONV - g^g convection heat 
loss, QRAD = radiative heat loss) , 
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heater temperature at that point), k is die thermal conductivity of the gas phase, t 

s 

and Nu is the Nusselt number for the convective heat transfer. Based on some heat 

(9) 

transfer calculations, ' a value of Nu = 6 has been used in ir.ost of the calculations; 
in essence, the heat loss is assumed to be ~ 6 times that by strai^t conduction through 
the gas. The convective heat is not actually ca.*ried away by the flowing gas; rather, 
there is a heat exchange between the ribbon and ti.o afterheater block at each point. 

The Nusselt number is clearly an adjustable parameter in the calculatic.>s; its 
appropriate value for a given growth system should be obtainable from experimental 
measures of the thermal profiles and comparison to theory. 

The boundary conditions to the numerical calculations are that the interface 
temperature T(0) = T^^, and that the temperature at the cold end of the ribbon 
T (L) = T = constant. Fixing the cold-end temperature allows one to obtain a con- 
vergent solution more rapidly; the effects of T on the vertical profile are readily 
examined in the model. The numerical solution uses the Gauss-Seidel technique 
with accelerated relaxation; a converged solution is found after approximately 100 
iterations. 

The output of the numerical calculations includes a graphical printout of the 
temperature profile and a prediction of the maximum growth rate for the given 
"experimental" conditions. The temperature values at the incremental points can 
be directly input into a numerical program vhich calculates the thermal stress in 
the ribbon. Some examples of the calculated profiles will be given in the following 
to show the effects of specific heat transport in the ribbon and the effectiveness of 
gas cooling of the ribbon surface. 

c. Effect of specific heat transport 

(6 7) 

A significant refinement to the earlier analysis' ’ ' which gave 
Eq, (1) has been the inclusion of the specific heat carried by the moving ribbon into 
the numerical calculation of the vertical temperature profile. The effect of this term 
on the thermal profile turns out to be quite important, as illustrated in Fig. 10. Fig. 

10a shows tiie temperature distribution in a stationary (V = 0) ribbon 4 cm long, 0.025 
cm (~ 10 mils) thick for radiation into a 0*K environment, whereas Fig. 10b is the steady- 


tThe temperature dependence of the gas thermal conductivity (e.g.. Ref. 8) has been 
account^ for in the calculations. 
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Fig. 10. Calculated vertical temperatu’^e profiles in a 0.025 cm thick, 4 cm long 
ribbon for (a) Vg = 0 and (b) nTjp = 0.117 cm/ sec. The convective heat 
transport is included in the calculations. Radiation into a 0°K environment 
is assumed; the tenperatura at the far end of the riobon is fixed to 
400PK. The arrow denotes :he position of the growth interface on the 
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state profile for a ribbon growing at 0.117 cm/'sec (~2.76 in./min). The latter rate 
corresponds to the condition where the gradient in the liquid is approximately zero, 
and is therefore the maximum ribbon growth rate for a 0.025 cm thick ribbon. Tlie 
pre /iou3 theory' ’ ' which neglects the convective heat transport predicts, on the. 
other hand, that the maxin",um growth rate is 0.146 cm/sec. Table I lists the values 
of the maximum growth rate as a function of ribbon thickness based on the earlier 
theo ’7 and on tlie current numerical calculations. The latter values are seen to be 
approximately 20% lower in each case. The same result is obtained for cases where 
the environment temperature is different from 0“K; the maximum growth rates based 
on the thermal profiles in the stationary ribbon <i.e., V = 0) are typically 10 to 20% 

o 

higlier than the rates obtained from the actual profiles in the moving ribbon. 


Table 1. Maximum Ribbon Growth Rate as a Function of Ribbon Thickness 


t 

(V ) , cm/sec 

g'max’ 

(V ) , cm/sec 

' g'max’ 

(cm) 

(using Eq. (16) in Ref. 7 ) 

(corrected value)' 

0.01 

0.231 

0.179 

0.02 

0.163 

0.130 

0.025 

0.146 

0.il7 

0.03 

0.133 

0.107 

0.04 

0.116 

0.094 

0.05 

0.103 

0.084 


♦The corrected value is obtained by numerical solution of the steady-state tempera- 
ture profile; the model includes the convective heat transport (i.e., specific heat) 
by the moving ribbon. Radiation into a 0°K environment is assumed. 


d. Effects of afterheater block and gas cooling 

The effectiveness of the heat exchange betweav the ribbon surface 
and the aftv^heater block can be best’’ Mustrated by somespecific examples. Figure 11a 
shows schematically a typical "hot" afterheater block configuration. The hot portion of 
the block (1 cm in lengdi) has a linear temperature gradient from 1300°K to 1000®K; 
the temperatures below and above the block are 1300*K and 500®K, respectively. The 

21 ‘ 


i 



0.025 cm 


0.05 cm 


0.05 cm 


t = 0.025 cm 


m I — 


,0.25 cm 


T „ , °K 
env 


Fig. 11. Typical examples of hot (a) and cold (b) afterheater blocks. The 

thermal environment of the ribbon is determined by the positioning of 
the block with respect to the ribbon and the growtn interface, and by 
the temperature values below, along and above block. 








block is located 0.05 cm (~20 mils) from tlie ribbon surface on either side and 
0.05 cm above the growth interface. The temperature profile in a 0.025 cm thick 
ribbon at zero gro>Mh rate is shown in Fig. 12a for the case where heat exchange 
between die ribbon and the afterheater is by radiation only (i.e., no gas cooling or 
heating). Figure 12b shows the profile for combined radiation and argon gas cooling. 
The tendency of the profile to be more linear within the block in the argon case is 
readily apparent from the figures. From the zero growth rate profiles in Fig. 12, 
one can obtain an estimate for the maximum growth rates in these systems (for a 
0.025 cm thick ribbon). The values are found to be 2.36 in./min and 3.08 in./min 
for the no-gas and argon cases, respectively; using helium as the heat exchange 
medium, one obtains a maximum growth rate of 5.79 in./min (as explained earlier, 
the actual maximum rates are about 10 to 20% lower than the above values as a 
result of the convective heat transport in the ribbon) . 

An example of a typical "cold" afterheater block is shown in Fig. 11b. Here 
the block is assumed to have a gradient of 240“K/cm over a 2.5 cm length. The 
block is located 0.25 cm above the growth interface and 0.1 cm from the ribbon sur- 
face. Temperature values below and above block are 1300“K and 300®K, respectively. 
The temperature profiles in a 0.025 cm thick ribbon are shown in Figs. 13a and 13b 
for no gas and helium, respectively. In the latter case, the profile in the ribbon is 
seen to attain the linear profile in the block over the last 2 cm of the block. The maxi- 
mum growth rates from these V = 0 profiles are estimated to be 2.70 in./min, 

3.20 in./min and 4.90 in./min for no gas, argon, and helium, respectively. 

The impoitant outcome of this numerical analysis is that we have a means 
of manipulating the thermal profile within the ribbon, which is Important insofar as 
the control of the thermal stresses. Clearly, there can be a direct correlation made 
b^ween these calculations and the design of the experimental hardware. The effects 
of heat exchange between the ribbon and the block via a gas phase were demonstrated 
by the above examples; means of controlling the profile by varying the gas compo- 
sition, for example, are also apparent. From the growth rates predicted for the 
afterheatea*s in Fig. 11, it is felt that the program objective of 3 in./min growth 
rate should be readily achieved using relatively simple cooling designs. The numeri- 
cal program allows one to consider trade-offs between the growtli rate and the ther- 
mal profile desired for low stress. Of course, the theoretical calculations must be 
reinforced with experimental measures of growth rates, thermal profiles and the 
boundary conditions which alter the thermal model. 
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Fig. 13. Temperature profile In a 0.025 cm thick ribbon at zero growth 
rate for the cold afterheater block in Fig. 11b: (a) No gas in 
system; (b) Helium gas cooling. The arrows denote position 
of the block relative to the profile. 
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e. Horizontal and transverse temperature profiles 

The analysis thus far has neglected gradients in temperature along 
tlie ribbon width (i.e., the horizontal profile) and through the ribbon thickness (i.e., 
the transverse profile) , Clearly, more heat is lost at the ribbon edges because 
of the greater surface to volume ratio. Also, since hetit is lost at the ribbon sur- 
face by radiation and convection, there must be a gradi ant through the ribbon thick- 
ness. An analysis of the three-dimensional temperature distribution in the ribbon 
was made on an earlier NSF program The calculations showed that the transverse 
gradients and the horizontal gradients at tli e ribbon edges were much smaller than 
the gradients in the vertical tempei r.ture profile at a given point. For example, 
the typical maximum temperature variatici. across the ribbon thickness is to 
5®; the thermal stress associated with this variation is much less than that which 
results from the non-linear vertical profile. In a later section, we will show how 
the horizontal temperature profile can be varied purposely along the entire width of 
the ribbon in order to control the thermal stresses, 

f. Linear vertical profile 

Finally, a brief discussion is given here on the prospects of 

obtaining a linear vertical profile in the ribbon vhich, according to the qualitative 

analysis in Fig. 8 is expected to lead to a condition of no thermal stress. We consider 

the element of volume between T. and Ti 4 j < T. in Fig. 9b. Even if the element were 

CONV^ RAn ^ 

perfectly insulated (i.e., Q. = Q.^ = 0), the gradient at T. would be less than 
at Ts because of the 1/T depaidence of the thermal conductivity. If additional heat is 
lost by the element (by radiation or convection), d T/dx will be even more positive. 
Therefore, in order to obtain a linear gradient, heat must actually be put into each 
element of volume (a small net heat input results from considering the specific heat 
transport in the ribbon). This, of course, would limit the rate of growth as the con- 
vective and radiative heat losses from the ribbon surface at the hi^ temperatures 
help to carry away the latent heat which originates at the interface. It is important to 
bear in mind in designing the cooling,/h eating environment that a growth rate of 3 in./ 
min for a 0.025 cm thick ribbon requires a vertical gradient of ~2400®K/cm in the 
ribbon at the growth interface. Clearly, a linear gradient of this magnitude cannot be 
maintained in practice over much of the ribbon length; basically, the quantity of heat 
associated witii this gradient cannot be transferred from the ribbon at the lower 
temperatures. 
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4, Thermal stresses in the growth of ribbon crystals 

This section summarizes the current status of the theoretical and ex- 
perimental efforts relating to the thermal stress problem in EFG silicon ribbons. An 
excellent description of the theoretical formalism of thermoelastic problems is 
given in Chapter 8 of Ref. 10; the efore, after presenting some gene“''.l remarks, 
we proceed to describe here the methods of solution and results only for the ribbon- 
shaped geometry of interest. First, we present an analytical solution for the ther- 
mal stresses in die ribbon, and discuss the limitations of this solution. We then 
describe a numerical method for calculating the stresses, and present the solutions 
for some simple temperature distributions. Procedures are described for deter- 
mining the plane and magnitude of the maximum shear stress at every point in the 
ribbon, and for calculating the resolved shear stress in the slip planes and slip 
directions of the dislocations for a given growth oriaitation. Experimental approaches 
to the thermal stress problem and some experimental results are also described. 

a. General remarks 

The basic assumptions in the tiieoretical formulation of thermoelas- 
tic problems (see Chapter 8 of Ref. 10) are: that the temperature can be determined 
independently of the deformations of the body, that the deformations are small, and 
that the material bdiaves elastically at all times. Also, it is assumed thcl ilie body is 
isotropic, and that the Young's modulus and thermal expansion coefficient are inde- 
paident of temperature (the latter assumptions are not necessary in the numerical 
calculation of the stresses). As described in an earlier section, thermal stresses 
arise as a result of non-uniform temperature distribution*? in the body (as well as 
stresses caused by external constraints). In the general three-dimensional problem, 
there are 15 unknowns; the three displacements (u, v, w), the six strain com*- n.ents 
(c .j), and the six stress components (ojj)* These variables are related by 15 equ - 
tions: three equilibrium equations (i = x, y, z) 


®”ix . °*’lv . 

3y 3z 

,th 


+ F^ = 0, 


(3) 


vdiere F. is thei compo' ent of the body forces; six stress-strain relations 
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and 


^xx = I [""xx ■ ^ ^zz^] 

^xy" (~^) “xy’ 

where E is Young's modulus, p is Poisson's ratio and a is the thermal expansion 
coefficient; and the six usual strain-displacement relations (e.g., p. 248 of Ref, 10). 

In addition, boundary conditions (surface tractions or displacements) are specified 
on the bounding surfaces. It should be noted from Eqs. (4) and (5) that the thermal 
expansion contributes only normal strains (i.e., no shearing strains). 

A principal prediction of thermoelasticity theory is that if the temperature 
varies linearly with a set of rectangular Cartesian space coordinates, then all the 
stress components are identically zero throughout the body (provided that surface 
tractions, body forces and displacement discontinuities are absent). Conversely, 
this is the only temperature distribution for vhich all stress components are 
identically zero. This concurs with the qualitative argument we presented in Fig. 8, 
The general, three-dimensional thermoelasticity problem cannot be solved 
analytically; various methods of solving for the strains, stresses or displacements 
are described in Ref. 10. However, the problem of calculating the stresses in the 
ribbon is considerably simplified by recognizing that the concept of plane stress is 
applicable in this case. Here we have the temperature T = T (x, y), and the thickness 
t in the z-direction small with respect to the other dimensions (see Fig, 14); therefore, 
it follows that 


*^zz “ *^xz 



( 6 ) 


The remaining stress components are obtained from the stress function 0 (x,y) 
which is a solution of the following differential equation: 


V4 0(x,y) 


3^0 . „ , 3^0 

^ 3x^ 3y^ 3j^ 


= “Ea 


W 



(7) 


2a 



thickness = t 


! 



Fig. 14. Definition of coordinate system used in the calculation 
of the thermal stress in the ribbon. The plane stress 
solution is applicable for t « C and t << L, where 
L is the ribbon length. 
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The boundary conditions for traction-free surfaces (which is the case of interest to 

00 

us) are given by 0 = ^ = O, where n denotes the normal direction to the bounding 
surface. The stresses are given by; 


and 


XX 


^xy 


0^0 

0y 




0^0 

9x 


0^0 

0\0y • 


( 8 ) 


In terms of the stress components, the boundary conditions are that - O 

2wv xy 

at X = O, and that a = o =Oaty = ±C, where C is one-half of the ribbon width. 

yy xy 

b. Analytic solution for the stress in the ribbon 

Boley and Weiner^^®^ present a solution for the stress function 9 
for the case of a thin rectangular beam (see Fig. 14) of width 2C and length L with 
C/L << 1, The temperature in the beam Is an arbitrary function of the coordinate y 
across the ribbon width, but it is a some\nAiat restricted function of x, as discussed in 
more detail below. The solution to Eq. (7) and traction-free surface boundary conditions 
is obtained as an infinite series 0 = 0^ + 0^ + 0g + ... ; successive terms of this 
series depend on successively higher derivatives of T with respect to x. For example, 
0. depends on 0^^^ ^^T/0x^^^ . Therefore, the series will converge rapidly if 

the temperature variation with x is sufficiently smooth, i.e., if it is expressible in 
a sufficiently rapidly convergent power series. The solutions for 0. up to i = 3 are 
given in Ref. 10 for the general case where T = T(x, y). 

When the temperature is a function of the vertical coordinate x only, the first 
three t£arms in the solution are given by 


and 


(v'-c2). 


dx 


(9) 


The leading terms for the stress components can be obtained from Eqs. (8) and (9): 
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•s 

aEC^ 

6 


(10) 

’“th 

oEC^ . 
6 


(11) 

th 

aEC^ 
24 ’ 


(12) 


where the subscript "th" is used to denote the tliermal component of the solution. 

From these equations, we see the qualitative d^endence of the stress components 
on the ribbon width (w = 2C) and on die derivatives of temperature. In the qualita- 
tive discussion of the origin of thermal stresses in Fig. 8, we were essentially con- 
sidering the <r component of the stress. 

(i) Corrections for non-zero tractions at interface 

One of the problems with the above solution is that the require- 
ment of zero traction at the end x = 0 (i.e., the solid-liquid interface) is not satisfied; 
rather the tractions are self- equilibrating there (i.e., no net forces or moments). 
According to Saint- Venant’s principle (e.g., p, 345 in Ref, 10), the above solution is 
therefore valid only at distances x >,2i,. However, most of the plastic deformation 
in the case of the ribbons occurs in the region x < 2c; hence we must add appropriate 
correction terms to the stresses in Eqs. (10) to (12) in order to account for the non- 
zero tractions at the interface. 

The problem of the stresses arising from the application of self- equilibrating 
tractions on one end of a long, rectangular strip has been treated in a series of 
papers by Horvay,^^^' In order to obtain zero traction at the interface in die case 
of the ribbon, we must apply at x = O the normal stress (cf. Eq. (11)) 

§} 

■ 


• (y)= - a^. 

6 ^1 


and the shear stress (cf. Eq. (10)) 


y. (y)= - d’T 

^ xy ^ V 6 


x = 0 




dx 


x= O 
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The stresses which result from the application of the normal stress in Eq. (13) are 

12 ) 

given by 




a EC* d*TI 

6 * 2 
dx 




a EC* d*T 


x/C 


dx' 


= 


aE d^T 

24 


dx 


. . sin (^gX/C) 

jcos OSgX/C) + ^ sin (^2x/C)j 


p 02X/C)J, 


(15) 


• e 


while the shear stress distribution in Eq. (14) leads to the stress components *** 


, , . oEC* . d*Tl 

‘<'xy> . 6- -T^ 


("icx' 

T 


a EC® d®T 
® dx® 


(cTyy) = - 


aEC* 


-L /i ^ ^~oioX./c • 

*C* -tj • ® 2 

= 0 ^ 
jcos (^gX/C) - sin 032^/<^)j 

•f (-<)• 

r=0 ' C*/ 

• e'“2^^^ . sin 

■ ^ /l if e" 

X = C 

/^2_^2\ 

[cos OjX/C) - > 


(16) 
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vtiere 012 = 2.0751 and = 1.1429. Similar correction terms for the higher order 
terms in the "thermal" solution in Eq. (9) have also been derived. 

The total value of each thermal stress component is thus the sum of three 
terms, viz.. 


ioii) -( 0 , 4 ) +( 0 . 4 ) 


ij 


total 


U 


th 


ij 


4 (.jj) 


(17) 


Therefore, each stress component is a more complex function of the ribbon width 
and of the derivatives than given in Eqs. (10) to (12). Ihe y-dependencies of the 
correction terms in Eqs. (15) and (16) are seen to be the same as those of the respec 
tive terms in Eqs. (B) to (12); the magnitude of the terms decrease exponentially 
along the vertical coordinate. An example is considered in the next section. 

(ii) Thermal stress for radiation into 0**K environment 

Consider the case of a ribbon which radiates into a 0*K en- 
vironment. The temperature distribution is given byt 


T(x) = 


(a X + b) 2 


where 


a = 


8 . (w Ul^Jl 

3 ^ k T 


m m 


and b = 1/T, 


m 


(18) 


We used the following values in the calculations: w = 2 cm, t = 0.025 cm, 

•5 2 4 7 

c - 5,67 XlO erg/cm (®K) sec, = 0.22 x 10 erg/cm ®K sec, c = 0.54 and 
Tj^ = 1685®K, The resulting temperature distribution is shown in Fig. 15; it is 
readily seen that plastic deformation in response to the thermal stresses can be 

expected only in the region x < 1 cm ^ere the temperature is in excess of -'600*0, 

. 

Figures 16 to 18 show the three stress components calculated for the tem- 
perature profile in Fig. 15. In each case, the contributions of the various terms in 
Eq. (17) to the total value of the stress component are shown. In the top pan of 
each figure, we show the x-dependence of the maximum of each stress component, i.e., 

<^vvl » latter case, the maximum 

”iy/C»l yy|y/C»0 ^y/C*0.6 


tThe various terms in these expressions are defined in a previous section. The 
ribbon is assumed to be stationary for this calculation. 
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/orE at y/c = 



Fig. 17. Thermal stress component oyy in the ribbon 
for radiation into a 0°K environment. 
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actuall}' occurs at y/C = 1/^7 ). The y-dependence of each stress component can be 
obtained by combining the top curves with the normalization curves in the lower 
parts of Figs. 16 to 18. 

Considering die stress component cr in Fig. 16, we see that the outer 
portion of the ribbon (i.e., j y/c|> 0.6) is in compression, whereas the inner por- 
rlon is in tension. It is likely that the plastic deformation during growth, in 
response to these thermal stresses exceeding the yield stress, accounts for the 
"lengthening" of tlie ribbon in the central region, as observed from the residual 
stress in the ribbon (cf. Section U.B.2.). 

(iii) Maximum and resolved shear stresses 

Ratiier than discussing the above results in terms of the 
individual stress components, it is more useful to evaluate tlie maximum shear stress 
at each point in the ribbon and to determine die plane on which this stress acts. 

As described in any standard text on elasticity, there exist two mutually per- 
pendicular sets of axes on which the shear stress vanishes; the normal stresses on 
these planes are die "principal stresses" and 02 which are given by 




rxx - <^5 


The principal axes are at angles 9 ^ and x-axis, where 


tan 9 


‘^xx"‘^l,2 


The shear stress is maximum on a plane midway between 9 ^ and ^Iso 

on a plane perpendicular to this plane. The maximum shear stress is given by 




Table II lists the values of t (in kg/mm ) as a function of position in the 
ribbon for the previous example of radiation into a 0®K environment; the plane on 
which this maximum shear stress acts varies from point to point in the table. t 


t Values of E = 1.94 x 10^ kg/mm^ and a = 5 x 10“® ®K ^ have been used in the cal- 
culations. 
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Table II. Maximum Shear Stress (in kg/mm ) as a Function of Position in the 
Ribbon, and Temperature and Yield Stress as a Function of Vertical 
Distance. 


X (cm) 

0 

y (cm) 

0.3 0.6 

1.0 

T 

(»K) 

(kg/mm^) 

0 

747 

618 

306 

0 

1685 

0.03 

0.2 

52 

83 

103 

45 

1359 

0.40 

0.4 

49 

53 

56 

81 

1170 

1.70 

0.6 

6.1 

9.8 

25 

89 

1043 

4.45 

0.8 

20 

22 

27 

81 

950 

9.2 

1 

31 

31 

31 

66 

878 

16.4 

2 

11.1 

10.9 

9.6 

10.0 

668 

83 

3 

1.03 

0.95 

0.72 

0.54 

560 

190 

4 

0.08 

0.06 

0.02 

0.17 ' 

492 

330 

5 

0.08 

0.06 

0.02 

0.17 

444 

470 


♦Yield stress values were estimated from J. R. Patel and A. R. Chaudhuri, J. Appl. 
Phys. 2788 (1963). Low temperature extrapolations probably exceed fracture 
stress. 

The y = 0 values correspond to the central axis along the ribbon, y = 1 cm is the rib- 
bon edge (a width of 2 cm was assumed). The table also lists the values of the tem- 
perature and of the yield stress of silicon at each temperature. It is seen that the 
thermal shear stress (if the material were to bdiave purely elastically) exceeds the 
yield stress in most of the region x ^ 1 cm, and hence plastic deformation will occur. 
The prediction of the residual stress in the grown ribbon is beyond the scope of the 
present analysis. Clearly, time dependent plasticity effects must also enter the 
problem; each portion of the ribbon at x = 0 must also experience the thermal stress 
at X > 0 as the ribbon cools down to room temperature. 

Since dislocation glide is likely to be the predominant mechanism of plastic 
deformation, one should also calculate the resolved shear stresses on the slip planes 
and along the slip directions of the dislocations for a given crystal growth direction. 
Table III lists the resolved shear stresses for a ribbon growth orientation of (111) 

[112]. For randomly polycr ystalline ribbons, such an analysis, of course, would not 

fl* 20) 

be very useful; recent evidence' ’ ' suggests, however, that the ribbons attain a 
{110} <112> growth orientation and a steady-state defect structure. 
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Table III. Resolved Shear Stresses for (111) [112] Ribbon Growth 
rit^ne L>ir action Resolved Shear Stress 


(Hi) 

<iio> 

4 

yrr 


(111) 

<011> 

4 

yio2 


(111) 

<101> 

4 

>/l02 

^^xx ^ °^xy^ 

(ill) 

<011> 

2 

s/V 


(ill) 

<ioi> 

2 

>/90 

{3c -(j* +20* ) 

XX ^yy xy 

(iii) 

<110> 

2 

3 / 2 ~ 

‘"xx "xy* 

(lii) 

<101> 

2 

z/V 

1 

Q 

(lii) 

<011> 

2 

>/90 

^xx ■ ^yy ■ 2 ‘'xy^ 

(lii) 

<110> 

2 

z^fY 

^*^xx " *^xy^ 


(iv) Limitations of the analytic solution 

The major problem with the analytic solutions described in 
this section is that convergence of the ' \ occurs only if T(x) can be expressed 
as a polynominal in x. In the derivation of Eqs. (10) to (12), which we used in 
subsequent analyses and examples, we have already neglected the higher order 
terms (l.e,, i > 3) In the solution. Thus we have neglected a d T/dx term in 
which would obtain from V*- in Eq, (9), whereas we included the same deriva- 
tlve in the calculation of a (from in Eq. (9)). It turns out that the tempera- 

yy ^ 

ture distribution (Eq. (18)) for radiation into 0“K environment actually leads to 
increasingly higher values of the derivatives; for example, the thermal contribu- 
tion to a from the i = 3 term is nearly an order of magnitude larger than the 
i = 2 term \diich we considered. Accordingly, the analytic approach outlined above 
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is only useful for temperature distributions vbich can be expressed in cubic (or less) 
power series in x. Otherwise, the solution gets even more complicated, as correc- 
tion terms must be added for the hi^er order terms in a and a in order to 

XX 

have zero resultant tractions at the interface. 

A further problem with the analytic solution is that it is not readily com- 
patible with the numerical calculation of the temperature profile in the ribbon 
described in a previous section. For these reasons, we have developed a numerical 
program for calculating the thermal stress in the ribbon; this will be described in 
detail in a later section. 

(v) Effect of horizontal temperature distribution 

The analytic solution^^^^ extends to the case where the tempera- 
ture is a function of the coordinate y along the ribbon width. Some preliminary calcula- 
tions werethereforemade to seethe effect of curved isotherms on the thermal stresses. 
The following form was assumed for the temperature distribution in the ribbon; 

T(x,y)=(l-,8y^).T(x). (22) 

Thus 0 > 0 corresponds to concave ("frowning”) and |3 < 0 to convex ("smiling") 
isotherms. The calculations showed that /3 > 0 leads to lowe r sti’esses, while 
8 < 0 leads to higher stresses than the 8=0 case. An approximate calculation 
also showed that the isotherm radius of curvature (with 8 > 0) must be of the order 
of the ribbon width in order to cancel the thermal stresses caused by the x-d^endence 
of the temperature. A simple qualitative argument (as in Fig. 8) readily shows that 
the frowning isotherms should indeed lead to a reducticr; in the stress component 
a caused by the vertical temperature profile. 

XX 

In view of these results, it is Interesting to speculate on how the growth of 
zero dislocation density crystals is achieved in the Czochralski technique. There 
the growth Interface is reported to attain a concave shape under the right growth 
conditions. It is conceivable that the combination of the radial temperature varia- 
tion (i.e., frowning isotherms) and the smoother axial profile in the crystals 
leads to a condition v^ere the thermal stresses during growth are negligible. 

There has been no general treatment of thermal stresses in Czochralski crystals 
caused by the combined effects of radial and axial temperature profiles. 



c. Numerical calculation of the stress in the ribbon 


S 


It is apparent from the above that it will not be possible to obtain 
an analytic expression for the thermal stress in the ribbon which is of such a general 
form that it could deal with the temperature distributions of practical interest. We 
have therefore developed a computer program for determining the stresses by the 
method of finite differences. 

At the point (I, j) in the ribbon (see Fig, 19), we can write Eq. (7) in the 
following finite difference form; 

20 [^> (i,j)|-8[0 (i,j + l) + ^ (i + 1, j) + 0 (i,j-l) (i-1, j)| 

+ 2 |0 (i + 1, j 4 1) 4 ^ (i + 1, j-1) 4 <P (i-1, j-1) 4 (i-1, j 4 1)] 

4 0 (i, j 4 2) 4 (i 4 2, j) 4 ^ (i, j-2) 4 <P (i-2, j) 

= -aE6^ [t (i, j 4 1) 4 T (i 4 1, j) 4 T (i, j-1) 

4 T (i-1, j) - 4 T (i, j)) . (23) 

Here we assumed that the nodal points are spaced equally in the x and y directions, 
i.e., that 6 x = 6 y = 6 , For example, we consider a ribbon 2 cm wide X 5 cm long, 
and set the increments 6 = 0.1 cm. Because of symm^ry about the x-axis, we need 
to consider only half of the ribbon (see dark solid lines in Fig. 19). Since we would 
like to use the general form of the nodal equation (Eq. 23)) for all nodal points, we 
need to "extend" our grid beyond the ribbon (solid lines) in Fig. 19, Let’s consider 
this extension by examining the boundary conditions. 

Along the y-axis (i = 2), we have o =0 and a - 0. Therefore, 

9 9 ^ 

a^?^/9y^ = 0, d<P/dy = A and ^ = Ay 4 B. Also, d^<P/dKdy = O and a0/3x = C. 

The constants A, B and C can now be chosen arbitrarily, as this will not affect the 
magnitudes of the stresses. Therefore, we set A = B = C = O, \y4iich leads to 
<P (2, j) = 0 and (1, j) = 9 (3, j) for all j. Similarly, we have 9 (i, 13) = O and 
9 (i, 14) = 9 (i, 12) for all i. We also know that 9 must be symmetric abcnit the 
x-axis; therefore, ^ (i, 1) = ^ (i, 5) and (i, 2) = (i, 4). 

At the far end of the ribbon (x = 5 cm), we suspect that 9 must be some 
smooth function of x. We therefore impose the boundary condition that a... ~ con- 

2 2 J. J. 2 

stant, i. e., d 9/9k = A. Therefore, 99/8x = A x 4 B and 0 = Ax 4Bx4C. We 
solve for the constants A, B and C by fitting the above quadratic equation to the 
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1 2 3 4 5 6 7 8 9 10 11 12 13 14 

(y = 0) j 

Fig. 19. Finite element grid for numerical calculation of the 
thermal stress in the ribbon. The dark IJnes denote 
one-half of the ribbon; J = 3 is the center axis, 
j = 13 Is the ribbon edge, and i = 2 is the growth 
interface. 
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points i = 48, 49 and 50 (for each given j), then we extend the solution to ^ (51, J) 
and ^ (52, j). We can now write the nodal equation (23) for all i = 3 to 50 and j = 3 
to 12. 

(i) Method of solution 

The solution of the nodal equations requires the specification 
of the temperatures T (i, j) and an initial guess solution of the stress function ^ (i, j). 
We can then substitute a new value of ^ (i, j) during each iteration, viz.. 


9 

new 


(i, 




(24) 


Here 9 ,, is the value of the stress function before the iteration, and 9 ^ is the 
old func 

value of ^ (i, j) obtained from the nodal equation (23). R is the relaxation constant; 

R > 1 implies over-relaxation, and generally more rapid convergence. In our case, 
we found that R = 1.95 resulted in the fastest convergence; for R > 2, the solution 
diverged. The operation indicated in Eq. (24) is carried out point by point (j = 3 to 
12) and row by row (i = 3 to 50) during each iteration. It was found that the value 
of ^ (i, j) remained constant to 5 significant figures after 5000 iterations. Ihe 
stresses at (i, j) are then calculated from finite difference formulas equivalent to 
Eq. (8). 


The form of the initial guess solution did not affect the results, except per- 
haps the speed of convergence during the first ~ 1000 Iterations. Therefore 
9 (i, j) = 0 everywhere Initially was used in most of the calculations. 


(ii) Examples of numerical calculation of stress 

The initial attempt to calculate stresses by the above numeri- 
cal procedure was for the case of a ribbon radiating into a 0®K environment. As pointed 
out above, the analytic calculation gave only approximate results in this case 
because of neglecting the higher order terms in the solution. A comparison of the 
numerical and analytical solutions for this example shows that the qualitative nature 
of the analytical solution (i.e„ the x-y spatial dependence and the relative magni- 
tudes of the stress components at a given point) is correct, but that the stresses 
shown in Figs. 16 to 18 are too large by about an order of magnitude. This does 
not imply that the numerical solution is necessarily correct, however. 
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In order to examine the accuracy of the numerical calculation, we need to 

consider a temperature distribution for vdiich the analytical procedure gives exact 

results. For simplicity, we assumed a quadratic form for T(x), i.e., T(x) = A + Bx 
2 

+ C X (e.g., A = 1685, B = -560 and C = 56 were used). A comparison of the thermal 
stress components obtained by the numerical and analytical techniques showed that 
the former were consistently with;:i < 10% of the correct analytical values. A 
sli^t problem with the numerical solution appears at the far end of the ribbon 
(x ^ 4 cm); here the numerical solution begins to diverge from the analytic solu- 
tion after many iterations. The problem undoubtedly lies in the far-end boundary 
condition used in the numerical calculations; this will require further investigation. 

In conclusion, it is apparent from these results that, vdiereas the analytic 
solution was useful to show the qualitative dependence of the thermal stress on the 
spatial variables and the ribbon width, the numerical approach is far more useful 
for the practical cases of interest to us. For example, the numerical program can 
determine the stresses for any given T (x,y) in the ribbon, such as may be found 
in a given growth setup. The accuracy of the solution (~10%) is clearly sufficient 
for our purposes; some effort must still be made to improve the far-end boundary 
condition. The compatibility of the numerical stress calculation program to the 
thermal profile program described in Section II.B.3,is readily apparent; the ultimate 
usefulness of these programs would be, of course, to determine the T (x,y) (and the 
appropriate thermal environment which leads to it) such that, on the time-scale of 
EFG ribbon growth, no plastic strain occurs at any point in the ribbon. Although the 
thermoelasticity calculations described here are not sufficient to predict the residual 
stress in the ribbons, they do show the extent of plastic deformation which occurs and 
the temperature and ribbon regions where plastic flow is most severe. The inclusion 
of time-dependent plastic flow effects into the model would, in principle, allow one to 
calculate the residual stress in the ribbon. At the present time, there is not sufficient 
data on plastic flow in silicon to permit such a calculation, 

d. Experimental approaches to thermal stress problems in EFG 
ribbons 

The foregoing theoretical analyses of thermal profiles and stresses 
in ribbon- shaped crystals have greatly influenced the experimental program, particularly 
the design of the experimental equipment used to control the thermal environment of the 
ribbon. The detailed results of the ejqjeriments were presented in earlier reports on 
this program; asummaryofthegrowthrunslsglvenelsewhereinthisreport. 
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(i) Control of the vertical temperature profile 

The vertical profile in the ribbon impacts both program objec- 
tives of fast and stress-free ribbon growth. For high crystal growth rates, we need to 
attain a hi^ initial vertical gradient (i.e., at the growth interface) in the ribbon. 

This requires rapid cooling of the ribbon surface in this region, which will then 

lead to a large value of the second derivative of the vertical profile. The possible 

2 2 

necessity for a trade-off between fast growtii (large dT/dx and d T/dx ) and low 

2 2 

thermal stresses (small d T/dx ) is therefore an immediate consideration. 

The control of the vertical profile has been implemented by various passive 
and active afterheater block designst (cf. Fig. 11). The extent of the hot (or cold) 
zone, the temperature profile along the block, and the location of the block with 
respect to the ribbon (i.e., distance from ribbon surface, and height above growth 
interface) have been the major variables in the experiments. Cooling of the 
blocks (by water flow) was also included in the designs. Typical growth rates vhich 
were achieved using "cold" afterheaier blocks (e.g.. Fig. 11b) were in excess of 
2 in./min and, in some cases, approached the program goal of 3 in./min. These 
rates are not any higher, however, than vhat can be achieved by simpl) shielding 
hot surfaces and allowir ; the ribbon to radiate into an essentially 0*K environment. 

In addition, all the ribbons grown from the cooled blocks were excessively stressed. 
The use of insulating blocks^^^^ has resulted in lower growth rates (~1.2 in./min to 
1.8 in./min) and still large residual stresses. Actively powered afterheater blocks 
have been designed and will be tested. To date, no afterheater block configuration 
has been used on this program which permitted growth rates In excess of ~1,2 ln./min 
and resulted in tolerable residual stress in the ribbon. Later in this section, we wlU 
describe a recent design strata \»hich is expected to lead to fast, stress-free ribbon 
growth; the design makes extensive use of the theoretical modelling in Section II.B.3. to 
obtain successive manipulation of the vertical temperature profile in the ribbon. 

(il) Control of the horizontal temperature profile 

The concqpt that concave (or frowning) isotherms in the ribbon 
should lead to a reduction in thermal stresses has been verified, to some extent. 


fA "passive" block receives its heat from the environmo’t (susceptor and crucible, 
dieheaterc, ribbon, etc.); an "active" block contains an independently powered 
heating element. 
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in some fast ribbon growth experiments on an earlier in-house program/ ’ It was 
found that tiie tendency of the ribbon to buckle (cf. Fig, 7a) could be correlated to the 
shape of the solid-liquid interface. The latter was varied by changing the radius of 
curvature of the die top (in the plane of the die). At steady-state growth, the inter- 
face was either smiling (in the case of .12 in. radius of curvature dies) or frowning 
(in the case of flat-top dies); an example of each is shown in Fig. 20. The buckling 
of the ribbons after growth occurred with only the curved dies, i.e„ with the smiling 
growth interface; the tendency to buckle can probably be associated with a larger 
residual stress in these ribbons. Of course, the shape of the interface is also that 
of the isotherms in the near vicinity, at leitst, of the interface. Hence the lower 
residual stress in the ribbons with the frowning isotherms provides support for the 
theoretically expected d^endence of the fhei'mal stress on the horizontal profile. 

Attempts to vary the horizontal temperature profile (and to attain frowning 

isoiherms) during the current ejqjerl mental program have been two-fold; neither 

approach has been particularly successful thus far in reducing the residual stress. 

First, the use of reverse- curved (l.e., frowning) dies has been considered; initial 
(14) 

experiments' with a 12 in. radius reverse-curved die encountered difficulties in 

attaining full-vndth (1 in.) growth, however. The residual stress measured in one 

(14) 

0,8 in. wide ribbon grown from this die was ~ 25, 000 psi; the small reverse cur- 
vature was clearly insufficient to affect the thermal stress in the ribbon (cf. Section 
Il.B.4.b) . A second approach has been to modify the cooling block to obtain preferential 
cooling near the edges of the ribbon. Althou^ the horizontal isotherms were clearly 
affected, as evidenced by the fracture morphology of the seed or ribbon orystals, the 
residual stresses were still excessive. Clearly, further analytical work is required 
to identify the proper horizontal profile which would offset the effects of the vertical 
profile on the thermal stress in the ribbon. 

(iii) Annealing of the residual stress 

An obvious solution to the residual stress problem would be to 
subject the ribbons to a process of annealing at a hi^ enough temperature such that 
stress relaxation, bj' the diffusion-controlled climb of dislocations, for example, can 
occur. Previous experlence^^®^ with moderately stressed (i.e., < 20,000 psi) ribbons 
has indicated that the residual stress could be reduced to a negligible level by annealing 
the ribbons at 850*C for ~25 minutes. Similar experiments with the highly stressed 
ribbons grown on the current program has thus far proved to be unsuccessful; the 
ribbons would invariably shatter on be^ng loaded into the annealing ftirnace. The 
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fracture of these ribbons is probably caused by the additional mechanical and ther- 
mal stresses during the loading operation. 

Assuming that the time-deoendent relaxation process, in EFG libbons con- 
taining dislocations, is likely to be the diffusion-controlled climb of dislocations, 
we can make an estimate of the annealing times, at various temperatures, which 

would be required to achieve complete stress relaxation. Using the experimental 

(11 181 

data at 850“C r^jorted above, and possible activation energies' ’ ^ of 2 eV, 4 eV 

or 5.1 eV for the relaxation process, we obtain the relaxation times given in Table IV. 
It is likely that *’he true activation energj' falls somevtiere between tlie limiting 
values we have assumed. It is readily seal from the table that stress relaxation 
occurs very rapidly at temperatures above ~ 1200®C and that it is very slow at 
temperatures below ~700“C. These findings have been instrumental in the develop- 
ment of a new strategy, described below, for reducing the residual stress in the 
ribbons during growth. 


Table IV . Stress Relaxation Times as a Function of Temperature for 
Various Activation Energies 


Stress Relaxation Time (min) 



Activation Energy 


T (“O 

2 eV 

4 eV 

5.1 eV 

600 

9400 

3.5 X 10^ 

9.2 X lO*^ 

700 

609 

1.5 X 10^ 

8.6 X 10^ 

800 

66 

172 

293 

850 

25 

25 

25 

900 

10.3 

4.3 

2.6 

1000 

2.2 

0.19 

0.05 

1100 

0.58 

0.013 

0.002 

1200 

0.18 

0.0013 

0.0001 

1300 

0.067 

0.0002 

0.00001 
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(Iv) A new strategy for minimizing residual stress in EFG 
ribbons* 

The thermal stresses that ar*ise as a result of non-linear tem- 
perature gradients in the ribbon during growth were discussed in detail earlier in this 
section. The tfieory described the forces that would tend to cause plastic deformation, 
and the dq)endence of the forces on the thermal profiles and on the ribbon geometry. If 
no plastic strain is permitted to occur during ribbon growth, then the above theory 
is immediately useful; it predicts that either the ribbon breaks in a region where 
the thermal stress is too high, or the ribbon emerges with no residual stress. On 
die other hand, if plastic deformation occurs, the calculation of the residual elastic 
strain in the grown ribbon requires a solution of the complete elastic-plastic problem; 
the latter involves the inclusion of the time-dqjendent plasticity effects into the 
theory. 

At hi^ temperatures, plastic deformation can occur by several distinct 
mechanisms; the important ones for this problem are the slip of dislocations and 
the climb of dislocations. 

For slip to occur, dislocations must be present, and for more than a very 
small strain, dislocation sources must operate. For stresses above some level 
which is conveniaitly (but incorrectly) referred to as the ‘'yield stress", deforma- 
tion by dislocation slip is essantially instantaneous; that is to say, on the time scale 
of EFG ribbon growth, the strain rate is fast enough for the stress not to rise above 
the yield stress. 

The process of dislocation climb allows plastic deformation to occur in 
response to any stress, however small. "Tie rate of strain for a given stress, 
however, decreases very rapidly witli decreasing temperature. For the present 
problem, this means that the stress in the ribbon does not remain at the yield stress, 
but is spontaneously relieved by the climb mechanism. This relief may be complete 
or partial depending on the time and temperature. 

Thus we see that residual elastic strain in the ribbon can be avoided if there 
is either complete relaxation of stress or no relaxation of stress (i.e., no plastic 
strain) during the whole of the change of temperature of each element of the ribbon 

t The collaboration in this work with B. Chalmers, A. R. Chaudhuri, M. C. Cretella, 

C. V. Hari Rao, D. N. o‘ewett, B. H. Mackintosh, A, I. Mlavsky and K. V. Ravi is 
gratefully acknowledged. The work described here was done on Mobil Tyco 
funded programs, and was not a part of this contract. 
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from rhe melting tempei ature to room temperature. We are therefore concerned with 
the minimum time required, at various temperatures, to allow complete relaxation, 
or the maximum time during which no relaxation occurs. From the relaxation times 
in Table IV, we see that stress relief in the ribbon would be effectively complete at 
any temperature above about 1200®C, and chat almost no stress relief would occur 
below about 700“C. 

The above considerations suggest that the vertical temperature profile in 
the ribbon, ideally, should consist of three zones. The first zone (from the melting 
temperature to ~1200*C) should consist of an initially large temperature gradient 
(~2400*C/cm for 3 in./min growth rate), and a very' rapid decrease of the gradient 
away from the growth interface. This is required in order to achieve the desired high 
growth rate. Plastic strain in this zone should occur fast enou^ for the stress to be 
completely relieved. Of concern in this temperature r«^ion is the minimization of the 
plastic strain induced imperfections. This can be achieved by minimizing the second 
and hi^er order derivatives of the vertical profile, or by count-iring the effects of the 
vertical profile with an appropriate horizontal profile. The second zone should be a 
region of constant temperature gradient in vhich the stress-free ribbon cools from 
~1200*C to below ~600®C; no stress or strain is introduced in the ribbon in this zone. 
The third zone should consist of a further decrease in gradient that is needed to reach 
zero gradient at room temperature. No plastic deformation can occur in this latter 
zone; therefore, if the ribbon is stress free when it enters it, it will again be stress 
free when it leaves this zone at room temperature. The crucial question in this third 
zone is the elastic (thermal) stress which is developed; if this stress does not frac- 
ture the ribbon, it will then disappear as the gradient becomes zero at room temperature. 

The previously developed numerical thermal profile calculations played 
a major role in studying the'feasibility of the above concqpi for the vertical profile 
in the ribbon, and in designing tlie hardware of the cooling/afterheating environ- 
ments used to achieve the desired profile. An example of a cooling/afterheating 
block is shown in Fig. 21a. The distance above the growth interface is 0.102 cm 
(40 mils), while the distance from the ribbon is 0.076 cm (30 mils), 'fhe thermal 
profile along the block consists of a 0.127 cm (50 mil) cooling zone at 473‘’K, a 
same thickness insulating zone at 973®K, and an actively powered afterheating zone 
with a linear gradient of 35.4®K /cm starting at 1373®K. The calculated temperature 
profile in the first 1 cm of the ribbon is shown in Fig. 21b, The ribbon is assumed 
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to be stationary (V = 0); thermal interaction between the ribbon and the block is by 
radiation and conduction through argon. A maximum growth rate of 3 in. /min (for a 
0.025 cm thick ribbon) can be predicted from tliis profile. Within the linear gradient 
zone of the afterheater, the temperature profile in the ribbon soon attains that in the 
afterheater block. 

The above concept for fast, stress-free ribbon growth has been tested and 

(19) 

verified on a MTSEC in-house program.' Growth rates in excess of 2 in./min 
and ribbon widths up to 2 in. have been achieved; the residual stress in the ribbons 
was negligibly low. The cartridge design, with the built-in cooling/afterheating 
blocks, used in these experiments has been transferred to the current JPL program 
and will be used in the initial wide growth runs (on Machine No. 1) . 

(v) Summar^' 

• Residual stress in EFG silicon ribbon crystals 
results from the plastic deformation the ribbons undergo at the high temperatures. 
The plastic flow is caused by the excessive thermal stresses in die crystal during 
growth. 


• Thermal stresses result from non-linear tempera- 
ture distributions in the ribbon. The principal contributor to the thermal stress is 
the vertical temperature profile; concave (or frowning) horizontal isotherms would 
tend to decrease the thermal stress. Transverse (or through thickness) gradients 
in the ribbon probably do not contribute to the thermal stress. 


• Thermoelasticity theory, applied to the ribbon 
geometry, can predict the magnitude of thermal stresses and the extent of plastic 
deformation at each point in the ribbon. The theory cannot predict the residual elas- 
tic strain and stress in the ribbon; die latter involves the inclusion of time-dependent 
plasticity effects into the theory. 


• Stress relaxation occurs rapidly in EFG silicon 
ribbon at temperatures in excess of 1200®C; almost no relaxation occurs below 
~700“C. 


• The effect of growth rate on residual stress is two- 
fold: it affects the vertical profile via the specific heat transport in the ribbon, and 
it influences the extent of stress relaxation which can occur at a given temperature 
(cf. Table IV). 
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• The thermal stress in the ribbon during growth, and 
hence the residual stress, increases with increasing ribbon width; this is in agree- 
ment with experimental observations. 

• The ribbon thickness affects primarily che vertical 
profile in the ribbon (dT/dx is proportional to t ' ) . 

• The possible effects of the thermal stress on the 
steady-state defect structure^^^^ and on the density of plastic strain induced imper- 
fections, and the effects of solutes (carbon, impurities, etc.) on the stress have not 
beat considered in d^ail. 


• The theoretical numerical modelling program for 
the vertical profile in the ribbon has been used extensively to design a cooling/after 
heating block for fast, stress-free ribbon growth. Initial experiments using this 
design have proved to be successful in that growth rates of > 2 in./min and ribbon 
widths up to 2 in, have been achieved with only negligible residual stress in the rib- 
bons. 


C. iTiin Ribbon Growth 

One of the goals of the program is the growth of ribbons of controlled uniform 

tiiickness down to 0.1 mm. The impetus for establishing such a goal is based upon the 

dual requirements of low material usage (low cost) and sufficient thickness of the 

fl4) 

silicon for adequate solar cell conversion efficiencies. This goal was reviewed' ’ 
in the li^t of the many known variables that influence ribbon growth, the reliability 
of crystal growth machine components, solar cell conversion efficiencies and the 
possible impact on yield and cost. The examination indicated that very thin ribbon 
growth is probably undesirable and the possibl e economic advantages of thin ribbons 
are significantly offs^ by other factors. The discussion of growth related factors is 
paraphrased below as an introduction to ‘he experimental effort and as a rationale for 
its method. Thus, ribbon thickness t is given by a known analytical expression^®’ 
of the form; 

t = f(tj. s. h^. 0^) , ( 25) 

where t^ is the thickness of the die top, s is the meniscus hel^t, h^^ is the effec- 
tive height of the growth interface above the liquid level in the crucible and 0^ is 
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the angle between the tangent to the meniscus and the growth axis. Among these 
parameters, is a material parameter being related to factors such as the local 
orientation of the crystal surface. Tlie other param^ers, t^, s and h^^^ that control 
ribbon tliickness are process or equipment related variables. Consequently, the 
sensitivity of ribbon thickness to these variables is of concern. 

Shape stability in growth can be associated with the magnitude of dt/ds. Low 

|dt/ds I will result in greater stability. The detailed analysis^ of Eq. ( 25) shows 

that the ribbon thickness is more stable as s or h decreases. 

eff 

A more explicit dependence of the uncertainty in t due to uncertainties in 
the other geometric or process related variables results from the equation: 


At = 


-^1 At + 

iL 

As + 

at 

Ah 4* 

at 

'd 

ds 

^eff 

"eff ^ 

0 


A9 


(26) 


At^, the change in die top dimension can arise as a result of machining tolerances, 
which would be more critical as t^ decreases for thinner ribbon growth. In addition, 
small distortions of the die top due either to thermal effects or due to the formation 
of SiC particles on the die top would also adversely affect thickness control of thin 
ribbons. In general, |9t/3t^j « 1 independent of other variables. Cons..<juently, 
any changes in t^ are immediately reflected as changes in t. 

The change. As, in the meniscus hei^t is associated with tenperature 

fluctuations both on a short term and long term basis. In addition, s is also a 

function of h Present range of control is As is ± 0.05 mm with long time prospects 
eff 

of ± 0.025 mm. 

h^.j^ changes continuously as the melt level in the crucible drops. With melt 
replenishmair, it should be possible to keep to within 1.25 mm. 9^ changes 
by about ± 1® and is largely related to local crystal orientation effects. 

The main contribution to At comes from the j 3t/0s | As term, with the next 
most important contribution coming from the 9t/3t^ term, which is ^ At^. As s 
increases (for given tj and h^^p and thus causes a decrease in t, the error At not 
only Increases due to the As term but also due to other terms. Increasing h^^^ also 
increases the As contribution to At, offsetting the decrease that occurs due to the 
other variables. The experimental approaches to thin ribbon growth are thus 
clearly guided by theory; the controllable empirical parameters are t^, s and h^^^. 
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t , was the first variable experimentally approached. Early in the program, 

( 21 ) 

one-piece dies were tested to veduce splaying of top edges of the dies' and there- 
by to achieve better dimensional cor.trol. The experiments were successful and a 
variety of one-piece dies were used during the effort to better understand and achieve 

'thin ribbon' growth. Open, closed, and plugged end one-piece dies were evaluated; 

(141 

the latter to effect a higher meniscus and improve edge stability. 

Characterization of the dies included examining the interaction of the dies with 

the melt, specidcally the formation and distribution of SiC. Silicon carbide particles 

/ 28 ^ 

in ribbons are loci cf striictuxal imperfections and hi^ impurity concentration, ' 
and are detrimental to solar coll performance. Early in tins program, a one-piece 
die was sectioned, polished, and examined by metallography. The die had beoi at 
melt temperature for about five hours and about 18 inches of ribbon was grown from 
it. A large number of SiC particles were found at the bottom outside of the die 
(Fig. 22a) . Very uniform carbide layers and infiltration of graphite by silicon were 
found in the capillary channels with very few acicular growths. Carbide layer thick- 
ness was ~0.02 mm and penetration depth was ~0.22 mm (Fig. 22b). However, 
acicular carbides were grown in the feed slot across the top of the die. 

Growth runs were made to evaluate dies with minimum t^. These dies had a 

d 

capillary slot width of only 0.013 cm and "knife-edges" for a total t^ of 0.013 cm. 
Vertical 0.09 cm capillary feed tubes are coupled to the slot by a 0.07 cm horizontal 
bore hole. The first die used filled only with great difficulty. Ribbon growth was 
very unstable - the resulting ribbon was ~14 mils thick. A post mortem showed a 
large effective die top thickness due to SiC growth all over the die (Fig. 23) . The 
subsequent run also used a thin t^ die. Filling was facilitated by partially plugging 
the 0.07 cm hole below the die slot. The seeds used were etched to ~ 5 mils thick- 
ness. The first ribbon slowly increased in thickness from ~0 mils to 8 mils over 
7 in. of growth. 

The subsequent six starts all bdiaved the same way. Growth rates were all 
around 0.6 in. /min. Faster rates were not achievable due to the freezing of silicon 
in the capillary. 

The above experience, corroborating many previous observations, indicates 

that the effective t^ will not be a precisely controllable parameter, growth rates froia 

a small t^ die will be low, and the thin ribbon goal must be achieved ;h rough control 
d 

of h „ and/or s. 
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Fig. 22. (a) SIC growth at the bottom of the die. Particles appear to 
have grown on the top of the residual melt and were loosely 
adhering to the graphite (die) wall. 

(b) Carbide layer and infiltrated zone in one capillary wall. 
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Two growth runs were made to evaluate the influence of h^^^ on ribbon thick- 
ness. However, meniscus height variations masked any h^^^ effects. Because fixing 
h and varying the meniscus hei^t is very straightforward, these experiments were 
next run. 

Three growth runs were made to specifically relate meniscus height to ribbon 
thickness, t^ and h^^j were recorded. (The die top dimensions of nine dies were mea- 
sured using an optical comparator. The mean total die top thickness (t^) was 14.1 
mils. The range was 3.2 mils and the standard deviation was 0,955 mils) . Die top 
thickness at the end of each run but prior to cool down was measured from calibrated 
Polaroid photomicrographs (Fig. 24a) . Dies of 3 radii were used to achieve three 
very different meniscus heists. The first run was mr de with a flat die; the second, 
with a 5" radius die; and the third, with an 8" radius die. Values for meniscus height, 
Table V, were also obtained from photomicrographs (Fig. 24b, c, d) . Values are 
corrected for the 14° viewing angle (actual meniscus height s = x.sec 14° - 1/2 tan 
14° (t^ - t) , where t is ribbon thickness and x is the apparent meniscus height) . 

^eff measured after the run was completed. 

Values for ribbon thickness wci-’e taken with a ball n.lcrometer measuring 
the center of the ribbon at the point indicated by the ribbon displacement indicator. 

The observed results are illustrated in Fig. 25 through 27. The experimental data 
demonstrates a variation of ribbon thickness with meniscus height, as expected from 
theory ( solid lines) . The quantitative disagreement between experiment and theory 
is probably caused by the difficulty in determining the edge-to-edge die dimension 
(t^) during growth and by the uncertainty in measuring the meniscus height from the 
photographs. A 2 mil change in die thickness, for example, can shift the theoretical 
curves to coincide with the experimental values. 

D. Characterization 

Silicon ribbons typically contain both crystallographic defects and discrete 
inclusions. The inclusions are clusters of SiC particles. The crystallographic 
defects are predominantly twins, intersecting twin boundaries, dislocations and 
low and angle grain boundaries. Experiments were carried out to detci mine 
if an optimum seed orientation exists which would give rise to single crystal ribbons, 
and to determine seed - ribbon orientation relationships. Experiments were also 
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(c) 


M) 


Fig. 24. Die top (.luring ribbon growth: (a) end view; (b, c, d) side 
view showing die, menUcus, and growing ribbon. 


Table V. Meniscus Height Versus Ribbon Thickness 

( mils) 


Ribbon No. 14-200* 


s 

t 

2.4 

14.5 

2.4 

14.5 

2.78 

15.0 

3.4 

14.5 

9.1 

10.5 

9.06 

9.5-10.0 

7.18 

11.5 

3.3 

13.5 

4.15 

13.0 

6.54 

10.5 

5.84 

12.25 

7.03 

10.25 


Ribbon No. 

14-201** 

s 

t 

8.35 

6.75 

6.71 

7.5 

7.45 

7.25 

8.18 

7.0 

7.00 

7.5 

6.13 

12 

8.78 

6.5 

7.54 

8 

6.71 

9 

8.81 

6 

7.41 

7 

7.44 

7.25 

4.10 

11 

5.06 

9.5 

6.71 

9.0 


Ribbon No. 

14-202*** 

s 

t 

8.99 

12.0 

11.77 

7.5 

11.32 

8.0 

6.90 

12.5 

7.56 

12.0 

6.8; 

13.0 

6.60 

12.0 

8.87 

11.0 

11.44 

12.5 

10.4 

12.0 

9.41 

11.5 


*hg^= 1.87"; die radius = *> = 20 mils. 

**h l.r>7"; die radius = 5", t . = 16.4 mils, 

erf d 

***h „= l.'H"; die radius = 8", t. = 16.5 mils, 
eff d 
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■e is the theoretical relationship; 






Ribbon thickness, mils 

Fig. 27. Ribbon thickness vs. meniscus height. Solid curve is the theoretical relationship; 
• denotes exper' mental results. 
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undertaken to examine the electrical activity of the various defects present in the 
as-groAvn ribbons. Spreading resistance measurements wore done on selected 
samples to assess the influence of the defect species on impurity distribution. 

1. Crystallinity and orientation studies 

Although single crystals have been grown by the EFG technique the 
silicon ribbons grown deviate from single crystallinity. However, the ribbons are 
not composed of variously oriented grains separated by hi^ angle grain boundaries. 
The sustained growth of long ribbons has been observed to result in an essentially 
uniform structure consisting of linear parallel boundaries largely parallel to the 
ribbon edges. 

The linear defects are postulated to be predominantly twins with some of 
the twin boundaries being associated with dislocations. The relatively hi^ densities 
of twin boundaries with the twin planes being normal to the ribbon surface should 
result in twin related regions between boundaries. However, X-ray data indicates 
that the ribbon surface has a predominance of a single orientation. Two possible 
mechanisms could account for this observation. If the twinned structure is com- 
posed of closely spaced twins or twin pairs with the crystal returning to its 
"untwinned" state immediately following the generation of a twin, the resulting 
structure would consist of a predominance of a single orientation with a small per- 
centage of the crystal possessing the twinned orientation. The relatively macro- 
scopic sampling procedure of the Laue back reflection technique or the channeling 
patterns obtained from an SEM mi^t not detect the twinned orientations between 
closely spaced twins. Some preliminary transmission electron microscopy studies 

have shown the presence of very closely spaced twins, with twin spacings as low as 
0 

250 A, in these ribbons. 

An alternative suggestion is that the linear boundaries are composed of 
stacking faults which would not give rise to twin reflections in X-ray micrographs. 

In the limit, a stacking fault is equivalent to a twin pair witli an extra plane of 
atoms (or two extra planes for the diamond cubic structure) between them (for the 
case of extrinsic faults). Figure 28 is a (110) projection of tlie diamond cubic struc- 
ture showing the crystallography of the twin pairs or an extrinsic stacking fault. 
Either of the above discussed models can iccount for the observations. However, 
attempts to more precisely assign one or 'he other of the mechanisms to explaining 
the phenomenon requires further investigation. 
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PLANES 


Fig. 28. {110} projection of the diamond cubic structure. The dotted 
lines r^resent twin planes. Notice the orientation relation- 
ships between the matrix and tlie twinned region. 
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Nucleation of the linear boundaries, whether they be twins or stacking faults, 
is very likely to be related to the presence of carbon in the material. The influence 
of carbon on the crystallographic quality of both melt grown and vapor deposited 
crystals of silicon has been well investigated.' ^ Carbon has been shown to 
generate growth twins, stacking faults and, at high levels of concentration, disloca- 
tions and high angle grain boundaries in silicon crystals and films. It has been 
suggested that twins observed in epitaxial silicon films are a result of the packing 
effects of two extra (COOl) SiC planes betw'^on the (111) planes of silicon. The 
incorporation of SiC in tlie form of short chains in the silicon lattice has been 
suggested as the nucleating mechanism for defects such as stacking faults. 

The closed circles in Fig. 28 suggest the atomic configuration that carbon atoms 
could assume in the ribbon resulting in the insertion of extra layers composed of 
carbon atoms wtiich then constitute the observed linear boundaries. Whether con- 
tinuous incorporation of carbon in the form of sheets of atomic dimensions in the 
growing crystal is required for the propagation of the observed defect structure is 
not known at present. It would indeed be sufficient for carbon or SiC to nucleate the 
twins or stacking faults with propagation occurring as a result of growth. 

Whereas the major defect types observed in EFG silicon ribbons are twins 
(or stacking faults), a small fraction of the structure is composed of linear arrays 
of dislocations often associated with the twin boundaries. Figures 29(a) and 29(b) 
show two examples of bands of dislocations. Frequently, deformation induced dislo- 
cations are also observed with dislocations piled up against twin boundaries. 

Two likely mechanisms can be invoked to explain the morphology of the dis- 
locations. Stress Induced or grown-in dislocations can be confined to regions 
between twins as shown in Fig. 29. The clustering of the dislocations away from 
the twin bands could be a result of image or surface forces of the twin boundaries 
exerting a stress on the dislocations resulting in the observed morphology. 

Stacking faults can be absorbed by twin boundaries by the displacement of 
twin planes with a displacement vector of g [112] during crystal growth.' This 
mechanism would result in stacking faults bound by partial dislocations of the type 
I {112} . A high density of overlapping faults in association with the twins would 
result in the observed high density of dislocations proximal to the twins. An 
analysis of the nature of the dislocations, i.e., vdiether the dislocations are of the 
type ^ {110} or partial dislocations bounding stacking faults, should shed more li^t 
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Fie. 29. Linear arravs of dislocations, often associated with twin 

boundaries. Arrays of dislocations between twin l>oundarles 
are stress Induced, as shown in (a). In addition, cluster^f 
dislocations, as seen in (a) and (b) are often observed. Th.e 
region in the immediate neighborhood of the twin is seen to 
be devoid of dislocations. 
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on the mechanisms. The attainment of a structure consisting essentially of parallel 
twins parallel to the ribbon edges suggests that following the establishment of steady- 
state growth conditions the defect stnicture of die ribbons essentially reaches a 
sieady-state equilibrium configuration. 

Orientation experiments involving variously ori mted seeds have been con- 
ducted with a view towards determining whether an "equilibrium defect" structure 
develops In EFG ribbons and if so what the effects of seed orientation would be on 
the establishment of such a structure. 

Table VI lists the seed orientations used as part of the investigation. Seeds 
v/ith (111) orientations were chosen because (111) planes are the pre/erreo twin 
planes in diamond cubic materials. I* has been reported that EFG silicon ribbons 
have a preferred surface orientation of (110).^^®^ (110) is a low index, hi^- 
symmetry plane. Hence, (110) seeds with different growth directions were chosen. 

A ribbon seed was also chosen for comparison with the (110) [211] seed. 

Nine ribbon growth runs comprised the growth orientation study - a total of 
32 growth attempts were made. The results are characterized below. For the orien- 
tation stVidy, the growth process was standardized as nearly as possible. The furnace 
arrangement included an active one-piece ribbon face heater (this, prior to the state- 
of-the-art conversion described earlier) as well as the crucible main iieater. A 
flat molybdenum heat shield sat on insulating stand-offs on the face heater. A two- 
piece graphite die with a 12 in. radius of curvature was used for all these orientation 
runs. Growth v/as from a 3 in. diameter quartz crucible contained in a graphite shell. 
Detailed growth t un data and results are presented in previous Technical Progress 
Reports of this program.^ 

The following observations were made during growth of these ribbons: the 
defect and surface structure of the ribbons are best where the ribbon is growing 
stably at full width, and at <'•0.8 in. /min; ( 110) [ 111] growth is very active right 
after seeding and the structure breaks down dramatically; (111) [110]; (111) [211]; 
(100) [110] liquid-solid interfaces are unstable and display facetting; (110) [211] 
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Table VI. Expected Ribbon Surface Orientations and Growth Directions 
after Twinning. 


Seed Orientation 

(110) [lii] 

Twinning Plane 

(111) 

(111) 

(ill) 

(lii) 

Orientation after Twinning 

(ii4) [151] 

(110) [iii] 

(114) [5il] 

(110) [Ii5} 


(111) 

(122) 

(4iii 


(lii) 

(122) 

[011] 

(100) [Oil] 

(iii) 

(122) 

[Oil] 


(ill) 

(122) 

[4ii] 


(111) 

(ii4) 

1712) 

(110) [112] 

(lii) 

(110) 

[ii2] 


(iii) 

(114) 

[172] 


(ill) 

(110) 

[552] 


(111) 

(ii4) 

[lio] 


(lii) 

(110) 

[ii4] 

(110) [lio] 

(iii) 

(ii4) 

[lio] 


(lii) 

(110) 

[ii4] 


(111) 

(iii) 

[lio] 


(lii) 

(151) 

[ii4] 

(111) [lio] 

(iii) 

(115) 

[lio] 


(lii) 

(511) 

[ii4] 


(111) 

(iii) 

[211] 


(iii) 

(151) 

[2i7] 

(111) [211] 

(iii) 

(115) 

[27ij 


(ill) 

(511) 

[255] 
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siuf^e crystal seeds most commonly lead to growth v'ith twins parallel to the growth 
direction and a smooth interface; the unstable 'orientation seeking* pe riod between 
initiation oi growth and steady -state growth is greatly reduced when seeded from 
EFG ribbons, effectively ( 110) [211] , containing twins parallel to each other and to 
the growth direction (the equilibrium defect structure described below) . 

A silicon crystal with a known surface orientation, growing along a given 
crystallographic direction shows first order twin traces on the ribbon surface such 
as those presented in Fig. 30; Table '/I lists the ribbon surface orientations and 
growth directions following twinning on the appropriate ( 111) planes for the selected 
seed orientations. Such twin traces have been seen in that part of the ribbon which 
is proximal to the seed. However, interactions b:tween primary twins leads to 
higher order twin generation, with a rapid breakdown of the simple matrix-twin 
relationships. The crystal achieves the ''equilibrium" state wherein a majority of 
the linear boundaries lie parallel to ea:;h other, and parallel to the direction of 
crystal growth. 

In order to establish the "equilibrium” orientation, y-ray Laue patterns were 
obtained at points on the ribbon, one foot and two feet from the seed. The orienta- 
tions were determined across the ribbon width at eight equidistant points at each 
of these regions. As an example. Figs. 31(a) and 31(b) show two X-ray patterns 
taken two feet from the seed in ribbons 14- 149-3 (seed orientation; (111) [ 110]) and 
J 4-154-3 (seed orientation; (111) [ 211]) From patterns such as those in Fig. 31, 
it was deduced that the predominant orientation of EFG silicon ribbon approximates 
to (110) [ 211] . Fig, 32(a) is a selected area electron channeling pattern (SAECP) 
from a region exhibiting the equilibrium structure ana is compared to a similar 
pattern from the (110) [ 211] seed in Fig. ;i2 (b) . Note the similarities in the 
diffraction patterns. This preferred orientation was observed in all the ribbons after 
about one foot of growth, irrespective of the seed orientation. Fig, 33 shows the 
orientation at the end of 2 feet of growth in sections of various seed orientations. 

Note the clustering of ribbon orientation near the (110) pole. The '.wsition at which 
the equilibrium structure is established was found to be dependent on the seed 
orientation. Figs. 34 and 35 show regions of ribbons near the seed for the (111) [110] 
and (111) [ 211] seed orientations, respectively. These show predomlnan'dy first 
order twin generation. However, the stable defect configuration is not achieved in 
these cases till much later in the growth process. In contrast, the stable orientation 
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( 110 ) 



[ 112 ] 






(111) ^ j 

[lio] 




( 110 ) 

[ill] 


Fig. 30. Schematic of principal tvdn systems in silicon relative to frequently used 
seed and growth orientations. (After Leipold, et al. Proc. IEEE Photovol- 
taic Spec. Conf. Rec, 11th, 1975) . 
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Fig, 33. Orientations of ribbons grown from the various seeds. The ribbc.i 
samples were two feet from .ced. Spots in stereographic projection 
correspond to the ribbon surface orientation. Notice the clustering 
near the (110) pole. 
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Fig, 34. First order twin generation in ribbon next to seed (orientation of 
seed: (111) [110] ). This seed orientation gives nse to inter- 
secting boundaries initially. However, twin reactions lead to ±e 
generation of parallel ivvlns and the equilibrium structure, as the 
crystal grows. Ar’*ow indicates seed-crystal interface. 
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Note diat there arc no parallel boundaries, parallel to tne crystal 
growth direction. Parallel boundaries form as crystal grows, and 
the equilibrium structure is eventually attained. N'^te also the 
generation of slip bands in both ribbon and seed. 
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forms immediately after seeding M^ien either the (110) [211] seed or a ribbon seed 
is used. Figures 36 and 37 show ribbon starts using such orientation seeds. 

In order to further examine the twinning orientations possible, a region close 
to the seed was chosen since the initial twin traces should correspond to first order 
twins making comparisons with theoretical predictions feasible. The seed orienta- 
tion was (111) [ 110] . Figure 38 is an example of such a region. It was found that 
the orientation of regions 1 and 2 was (111) , while region 3 had a (511) orientation. 
Thus region 3 is a first order twin obtained after twinning on the (111) planes. 
Regions 1 and 2 have the same orientation. This can be due to ( i) both regions 
being related in a matrix-twin relation, or ( ii) the linear boundary separating the 
two areas may be composed of multiple twins. A hi^ magnification photograph of 
the twinned region does not clearly reveal the multiple-twin nature of the boundary. 

In addition to twins and dislocations, discrete particles and clusters of 
particles of SiC are also observed in ribbons grown with the use of carbon d’es.^ 
These particles are postulated to be formed as a result of the rejection of super- 
saturated carbon by the solidifying liquid, the attendant precipitation of the carbon 
in the form of SiC on the die top and the subsequent incorporation of the particles 
in the growing ribbon as a result of the particles reaching the solid-liquid interface 
and being extracted by the solid. Although by employing a high meniscus, the 
incorporation of SiC particles can be minimized in the initial stages of growth, with 
sustained growth the particle dimensions eventually approach the dimensions of 
the meniscus and particle incorporation is resumed. The influence of the particles 
on the electrical properties of the ribbons has been investigated and it has been 
determined that the electrical activity of the particles is largely a function of die 
presence of lifetime reducing impurities around them. 

2. Electrical charaaeristirs 

With the relatively hi^ density of defects and the presence of lifetime 
reducing impurities the electrical characteristics of EFG ribbons are not comparable 
to conventional Czochralski crystals. Typical minority carrier diffusion lengths in 
these crystals range from 20 pM to 70 pM with typical Czochralski values being 
generally in e'xess of 100 pM. The importance of these numbers lies not in their 
lower magnitudes but in their relative non uniformity. The non uniformity of diffusion 
length will Impact the yield of the final product. In order to assess the contribution 
of the various defect types to the reduction in the diffusion length, a variety of 
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Fig. 36, Generation of parallel twin boundaries next to the (110) [211] seed. 
Arrows indicate seed-crystal interface. The (110) [211] seed 
orientation gives rise to pai'allel twins immediately after seeding. 
The few intersecting boundaries seen transform, as a result of 
twin interactions, into the equilibrium structure as the ribbon 
grows. 
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Fie. 37. Generation of parallel twin boundaries next to the ribbon seed. 

Arrows indicate the seed -crystal interface. As can be seen^fte 
linear boundaries in the seed propa^te into the growing ribbon. 
Observe some of the inter secring boundaries generating parallel 
twins. The structure very quickly attains the equilibrium orien- 
tation this case. 
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Fig. 38. Region of ribbon near the seed. Seed orientation was (111) [110]. 
Primary twinning on (111) planes can be seen. Regions 1 and 2 
have the same orientation, i.e., (Ill), while region 3 has the (511) 
orientation. These orientations conform to those expected for the 
given seed orientation. 
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experimental techniques have to be utilized. Details of ntany of tlie measurements 

( 29-32i 

made have been presented elsewhere.' In this section a brief discussion of 

tlie electrical activity of the linear boundaries and a disc ission of solar cell 
characteristics will be presented. 

3. Electrical activity of defects 

The influaice of the linear boundaries on electrical parameters and 
solar cell characteristics has been investigated to determine the degree to which 
the defects adversely affect die conversion efficiencies of solar cells.^^*^^^ A 
very effective method for the investigation of die electrical activity of crystallo- 
graphic defects in semiconductors is the use of the scanning electron microscope 
operated in the electron beam induced current (EBIC) mode.^^^^ This technique 
has been extensively utilized for die investigation of EFG ribbons.^^®* Figures 
39(a) and (b) show examples whereby the EBIC technique permits die imaging of 
electrically active defect boundaries. In Fig. 39(a) are shown secondary electron 
and EBIC images of a region of a solar cell fabricated in a silicon ribbon. The 
secondary electron display shows very litde structure on the ribbon surface 
whereas the EBIC image shows a number of dark line defects which correspond to 

linear boundaries which are electrically active, the dark contrast (or dark -white 
( 35) 

contrast) ' being a result of local recombination of Iieam generated carriers. 
Figure 39(b) shows a magnified EBIC image of another solar cell vha:e the two 
types of defect boundaries discussed in the previous section are identifiable. Linear 
twin boundaries displaying sharp contrast and the arrayed bands of dislocations 
which also function as regions of strong recombination are observed. The sharply 
delineated twin bands are also seen to be associated with dislocations. This fact 
is more clearly evideiit in Fig. 40. Figure 40(a) is an EfflC image of a region of a 
solar cell. A single linear boundary is observed to function as a region of strong 
recombination. Figure 40(b) is a secondary electron display of the preferentially 
etched surface of the same region as in ( a) . The arrows identify the electrically 
active boundary in the micrographs. In addition to the single electrically active 
botuidary, a hi^ density of electrically inactive boundaries is observed In (b) with 
none of them functioning as recombination centers as evidencei by their non- 
detectability with the EBIC mode of observation. Close examination of Fig. 40(b) 
permits the distinction between electrically active and electricsal’y inactive 
boundaries to be made. The electrically active boundary Is found to be composed of 
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(b) 



Fie. 39. (a) Secondary electron and EBIC images of a solar cell. Note 
the lack of structure in the secondary electron image ( picture 
on the left). The same area, in the EBIC mode, reveals 
marked contrast associated with regions of carrier recombi- 
nation, In ( b) are seen two types of defects that lead to strong 
recombination: twin bands with associated dislocations, and 
dislocation arrays. 
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Fig. 40. EBIC (a) ami srcomlary i-'k\'tmn (b1 images of a region of a 
solar cell. Only a single liiiL'af bc^uiit ai'v, as shown be arrov», 
is electrically active, (kliei lineat bouiKlaries are electrically 
inactive. A higlier magnification image of tlie electrically 
active boundary in (.:! shows that this boundary has dislocations 
associated with it. 
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dislocations possibly in association with twins, whereas no dislocations are 
observed in the vicinity or at the other linear boundaries. Figure 40(c) is a lugliei 
magnification secondary electron image of the electrically active boundri’ y showing 
the high density of etch pits due to dislocations. Both arrayed dislocations as well 
as discrete randomly dispersed dislocations have been observed to reduce the 
beam generated current in these crystals whereas the twins are largely 
electrically inactive. 

The strong recombination contrast observed at the defect boundaries is 

relab d to defect impurity association at these regions. Although clean, undecorated 

dislocations have been postulated to function as donors and acceptors in silicon 

it is very unlikely the., significant recombination can occur in the absence of an 

impurity at the dislocations. It has been documented that strong EBIC contrast is 

6s5“37) 

observed only when the defects are associated with impurities.' Consequently 

the dislocation arrays which are the primary electrically active defect species in 
the ribbons are undoubtedly decorated with a deep level impurity (or impurities) 
rendering them electrically active. 

With decorated dislocations and dislocation arrays being implicated as the 
principal lifetime reducing defect types in EFG ribbons, the avenues of approach to 
be followed for quality improvements cuc discernable. An overall reduction in the 
undesirable impurity content in these crystals should result in significant im- 
provement. In addition to reducing the impurity content of graphite dies the 
approach is largely one of increasing component and machine purity. A reduction 
in the dislocation density with particular reference to arrayed dislocations should 
also improve electrical quality. The methods to be adopted to reduce the dis- 
location density would largely be determined by the mechanisms responsible for 
the introduction of dislocations in these crystals. If dislocation generation is 
related to thermal stress in the crystals during growth, the establishment of 
appropriate thermal environments can reduce the dislocation density. On the other 
hand if dislocation introduction is a result of the presence of carbon in the 
crystals with perhaps tlie dislocations largely being partial dislocations in 
association with stacking faults the solution to the problem becomes much less 
tractable. Evidence to date, although insufficient, suggests that a significant 
fraction of the dislocations are stress induced. 
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An unanswered question is the effect of carbon in terms of directly influenci-'g 
the electrical behavior of EFG ribbons. Although carbon, unlike oxygen, is a neutral 
impurity in silicon, the high growth rates characteristic of the process may 
result in the occupancy of carbon atoms in other than substitutional sites, or the asso- 
ciation of carbon with oxygen, vacancies or other impurities resulting in electrically 
active complexes. Further work is required to assess the importance of such effects 
if they exist. 

4. Impurity distribution 

The large surface to volume ratio of the ribbons and the high vertical 
temperature gradients present during growth have an effect on tlie distribution of the 
dopant and impurities through the thickness of the ribbons. In addition, the position 
of the feeds slots in the die \^ich convey the liquid to the die top affect the distribu- 
tion of impruities and the dopant across the width of the ribbons.^ For steady- 
state growth with a quiescent liquid and sufficient temperature gradient in the liquid 
to a''Oid constitutional supercooling, each solute will have a concentration profile in 
the liquid from a maximum of C^/k^ at the Interface, decreasing exponentially to C^, 
the concentration in the melt. The concentration decreases to 1/e of its maximum 
value at a distance D/Vg from the interface, where D is the diffusion coefficient of 
the solute in the liquid and V is the growth rate. The question to be resolved is 

o 

whether the lateral flow of melt in tho growth film causes any significant lateral 
movement of the solute, either impurity or dopant, resulting in concentration 
gradients of the solute across the ribbon width or throu^ the thickness. Examina- 
tion of this effect for various solutes and various die geometries has shown that a 
heterogeneous distribution of impurities is possible in EFG ribbons as a result of 
liquid flow effects.^® ^ 

In addition to liquid flow effects, solid state diffusion of Impurities to surfaces 
(both external and inter/.al as represented by defects) is possible. By employing the 
surface photovoitage technique cf determining the minority carrier diffusion length 
in EFG ribbon solar cells, it has been determined that the minority carrier lifetime 
is generally different in the near surface regions of the material as comparexJ to 
that In the bulk.^ This suggests that an impurity gradient through the mickness 
can exist in these crystals. 

In order to examine the variations in concentration of the impurities along 
the width and through the thickness of ribbons spreading resistance measurements 
were undertaken. 
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A bevel was grounj and the surface polirhed using standard polishing tech- 
niques. The bevel angle used was 2“ 52', which leads to an apparent magnitication on 
tlic bevelled surface of 20. nic length of the bevelled region was typically 5 uO 6 mm. 
If the thickness of the ribbon is 0.2 to 0 25 mm (0.008 to 0.010 inches) , this would 
correspond to approximately 0.2 mm ( 0.008 inches) of the thickness of the ir.aterial 
removed. Two types of scans were made: 

1. along the ribbon width (approximately 2.54 cm), 

2. across the thickness and along the bevel, with L, C, and R 
referring to left, center and right scans respectively. 

The scheme followed is shown in Fig. 41. Probe steps of 10, 25 and 50 pm were 
employed; both single- and double-probe scans v, ere used. Carrier concentration 
was plotted, from Irvin’s curves,' ^ as a functioi of distance traversed. 

Figure 42 shows the carrier concentration in a ribbon grown from unbaked 

graphite dies. The ribbon was doped 'with boron to yield a resistivity of 1 n -cm, 

IB 16 

vdiich corresponds to carrier concentx*ations of 2 x 10 to 3 x 10 atoms/cc. The 
spreading resistance data shows the carrier concentration varying along the ribbon 
width by as much as two orders of magnitude. This is to be expected in silicon 
grown from impure dies, where a preferential segregation of impurities to structural 
defects may occur. Use of purified graphite components, reduces the fluctuations 
in resistivity (or, carrier concentration), as shown in Fig. 43. Small fluctuations 
in carrier concentration can be decreased by appropriate higli temperature anneal 
or gettering during solar cell fabrication. Fig. 44(b) shows the profile in a sample 
from the same ribbon as in Fig. 44(a), after hi^ temperature (~950®C) junction 
fabrication leading to tlie formation of a (phosphorus) layer. The junction depth 
was approximately 0.5 pm observed throu^ the tnickness of the material. 

In order to examine the factors responsible for variations of resistivity in 
the material, samples were examined by optical microscopy after final polishing 
and after a preferential silicon etch. The region marked on the right in Fig. 42 
was re-examined by using 10 pW probe steps (see Fig. 45a) ; the polished surface 
was also observed with an optical microscope. The surface relief in the optical 
micrograph of Fig. 45(b) occurs due to the variation in chemical polishing charac- 
teristics with resistivity. The presence of a grain boundary leads to an increase in 
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Fig, 41. Scheme followed for spreadii> resistance measurements. L, 
C and R refer to traces down me bevelled surface, i.e., into 
the sample bulk. 
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F*g. 42. Carrier concentration along the ribbon width. Ribbon was grown using unbaked 
graphite components. Probe steps are 50 pM apart. Note the large fluctuations 
in carrier concentration. 
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Fig. 44. Spreading resistance profile of as -grow’’, r'bbon. (a) shows 
the variation in carrier concentration, at three places, 
along the thickness of ribbon. Solar cell fabrication reduces 
the bulk concentration gradient, as seen in (b). 
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Fla. 45. (a) Magnification of spreading resistance plot of right hand 
region of Fig. 42. Probe steps are 10 )uM apart, (b) shows 
surface relief on as -polished ribbon surface, showing 
rfesence of a grain boundary, (c) Is the same area after 
preferential etching. Notice that not all grain boundaries 
cause a change In carrier concentration. 
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carrier concentration. Fig. 45(c) shows an optical micrograph of the same area after 
preferential etching. Note the presence of grain boundaries and dislocations, which 
do not seem to change the resistivity rtarkedly. Thus, not all grain boundaries act 
as gettering sites for impurities. However, grain boundaries tiiat do not change the 
resistivity appreciably (Fig. 46c) may act as effective recombination sites as evi- 
denced in Fig. 46(b), ’w^iichisa Schottky barrier EBIC scanofthe samearea, preferen- 
tially etched, as in Fig. 46(a) . 

It was observed that the carrier concentration both increased and decreased 
near grain boundaries. It is possible that such variations may be due to the differences 
in the nature of the grain boundaries. Tlie regions of a ribbon with hi^ dislocation 
densities showed a general reduction of the carrier concentration. 

E. Problems 
None. 

F. Plans 

Technical milestones for the last six months of this program include achieving 
growth rates of 7,5 cm/min (March 31, 1977) , growth of 7.5 cm wide ribbon (Feb, 28, 
1977), growth of 0.1 mm thick ribbon (Nov. 30, 1976) (Table VII). Emphasis on 
achieving the last of these goals is reduced. The programs to achieve high growth 
rates and wide ribbon described in the body of this report will continue on Machines 
JPL No. 2 and JPL No. 1 respectively. Concurrent studies of ribbon stress, solar 
cell characteristics, feasibility of melt replenishment and system purity will supple- 
ment the growth efforts. 


93 




Fic. 46. Comparison between defect morphology (a) Schottky barrier 

EBIC scan ( b) and spreading resistance ( c) . Crystallographic 
defects that act as effective recombination centers do not 
necessarily change carrier concentration* 
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APPENDICES 


1. An Updated Program Plan 

Table VII, unchanged. 

2. Man Hours and Costs 

Previous man hours were 3976 and cost plus fixe! fee was $142, 659. 

Man hours for July plus August were 1269 and cost plus fee was $48,587. Man 
hours and cost for September are estimated to be 659 and $25j 645 respectively. 
Therefore, third quarter man hours are estimated to be 1928 and cost plus fee 
to be $74, 227. Cumulative three quarter man hours and cost plus fixed fee are 
estimated to be 5904 and $216,886, respectively, 

3. Engineering Drawings and Sketches Generated during the Reporting 

Period 

Figures 47-53. 

4. Summary of Characterization Data Generated during the Reporting 

Period 

Spreading resistance data on ribbons grown from unpurified dies shows a 
large degree of scatter, whereas the use of purified dies results in a more uniform 
distribution of resistivity in the ribbons. Correlations have also been obtained among 
spreading resistance measurements, EBIC data and optical microscopy. 

5. Action Items Required by JPL 

None 

6. New Technology 

No reportable items of new technology. 

No cumulative new technology. 
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Fig. 47. Pedestal Base. 
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